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ION RESONANCE IN A MULTICOMPONENT PLASMA 


Ss. J. Buchsbaum 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received October 19, 1960) 


We have performed experiments designed to each ion species. At high plasma densities, curve 
study the phenomena associated with ion cyclo- B tends to the so-called electron-ion hybrid fre- 
tron resonance in a multicomponent plasma whose quency. This limit was studied by Auer et al.” 
charged particle density is sufficiently high that and in great detail by Korper,® who showed that 
interaction between the various plasma compo- for a plasma with a single ion species it is given 


nents takes place. We present here a prelimi- 
nary account of the experimental results. 

The transverse oscillations of a uniform, cold 
plasma in an axial, uniform magnetic field have 
been studied by a number of authors.'’~* In Fig. 1 = 
are shown the resonant frequencies of a hydrogen 
plasma which, for the sake of an example, is 1 
assumed to contain the ions H,*, H,*, and H,* in 
relative concentrations of 5%, 50%, and 45%, re- 
spectively. These frequencies were calculated** 
assuming no collisions by finding the roots of the 
equation ky(w)=, where ky is the propagation 
constant of an extraordinary wave propagating at 
right angles to the magnetic field. Equations (2) 
and (7) of reference 4 are examples of the results 
of such calculations. The curve marked A in 
Fig. 1 is given by 
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ELECTRON-ION 
HYBRID 


ww tw,, (1) 
where w p=eB/m is the electron cyclotron fre- 
quency and w» = (ne?/me,) is the electron plasma 
frequency. This resonance (included for com- 
pleteness) is of no concern in this experiment as 
it occurs at microwave frequencies and depends _ 
only upon the motions of electrons. Curves B, 10" = = 
C, and D are associated with ion motions and are NORMALIZED PLASMA DENSITY, @p?/ay? 
of interest here. At low plasma densities (“low” 
to be defined shortly), where sufficient space- plasma as a function of plasma density and magnetic 
charge field to couple the motions of the various field. The following ions are assumed present in the 
plasma components cannot be induced in the plas- indicated concentration: H*(5%), H,*(50%), and 
ma, resonance is at the cyclotron frequencies of H;* (45%). 


NORMALIZED RESONANT FREQUENCY, (@/Gp)? 
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FIG. 1. Calculated resonant frequency of a hydrogen 
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by 
w=w b” bi? (2) 

where w5; =eB/M; is the ion cyclotron frequency. 
For a many-component plasma w ; must be re- 
placed by DiX i bis where x; is the relative con- 
centration of the ith ion species. Curves C and D 
tend to the so-called ion-ion hybrids. It was 
shown in reference 4 in a three-component plasma 
that at these resonances the various ion species 
oscillate out of phase with each other while the 
electrons remain relatively motionless. The ion- 
ion hybrid frequencies are then sensibly inde- 
pendent of the electron cyclotron frequency and 
depend only on ion cyclotron frequencies and on 
the relative concentrations of the different ion 
species. For a plasma with more than two ions 
this dependence is sufficiently involved that the 
discussion of it must be postponed to a more de- 
tailed publication.® Here it suffices to say that 
the solid curves in Fig. 1 never cross the hori- 
zontal dashed ones, and that the larger the rela- 
tive concentration of a particular ion species the 
closer does its hybrid frequency approach the 
cyclotron frequency of the next lighter ion. Of 
importance in these experiments are the “knees” 
in curves B, C, and D at low plasma densities. 
They mark the transition from a single-particle 
behavior to a collective one; “pure” ion cyclotron 
resonance occurs in the region to the left of the 
knee. The density at which a knee occurs depends 
on concentrations of all ions,® but pure ion cyclo- 
tron resonance generally obtains provided 

Wy << W554) (3) 
where w ; is the cyclotron frequency of the ion 
in question. Ion-ion hybrids are attained at den- 
sities such that the inequality (3) is reversed, 
while for the electron ion hybrid one must have 
Wp? >> ws”. 

The aim of the experiment was to observe 
“pure” ion cyclotron resonance and to study the 
behavior of the resonance with increasing plasma 
densities. The experiment was performed on the 
positive column of an arc discharge with a therm- 
ionically heated oxide-coated cathode. The tube 
was connected to a vacuum system, which could 
be baked out and which could be filled with gas to 
any desired pressure. The gas used was “spectro- 
scopically pure” hydrogen, helium, or neon pro- 
vided by the Linde Company. Some, but not great, 
attention was paid to maintaining the gas purity. 
The positive column was in an axial magnetic 
field provided by a Bitter-type solenoid. The 


496 


large magnetic field (up to 85 kgauss) enables 
one to perform the experiment at elevated fre- 
quencies (6.2, 10.8, 18.18, and 27.3 Mc/sec were 


used) with a corresponding decrease in linewidths, 


Also, as (3) reduces to n<<[6x10"(B kgauss)?/ 
(M; atomic units)] cm~*, “pure” cyclotron reso- 
nance can be observed at densities which can be 
obtained in active gas discharges. The central 
portion of the positive column was surrounded by 
a resonant structure which, in the absence of a 
plasma, produced a “twisted” electric potential 
given by 


g(r, 0,2)=g) 7) cos(m6+k 2) coswt; m=1,2 


(4) 
where J,, is the modified Bessel function. The 
presence of the plasma modifies the argument of 
Im Which then exhibits the resonances shown in 
Fig. 1. The twist in the electric potential makes 
it possible for the rf field to penetrate the plasma 
more completely,® and is given by the axial wave 
number kz (kz >>w/c). The static magnetic field 
was uniform to within 4% over the interaction 
region. The absorption and dispersion were 
measured at a constant tube current and a con- 
stant gas pressure as a function of the magnetic 
field by measuring on an rf bridge the changes in 
the effective admittance of the resonance struc- 
ture caused by the plasma. This corresponds to 
motion along the diagonal lines in Fig. 1. No sub- 
stantial change in the absorption lines was found 
experimentally when the dipole field [m =1 in 
Eq. (4)] was replaced by the quadrupole field 
(m =2). 

Figure 2 is an example of “pure” or almost 
pure ion cyclotron absorption lines in helium with 
a hydrogen impurity. In Fig. 2(a) the H* and H,* 
lines are at their cyclotron fields while the H,* 
line is slightly shifted toward lower magnetic 
field, i.e., toward its hybrid field. There is also 
evidence for a shifted He* line. Figure 2(b) is a 
dispersion curve. 

In Fig. 3 we plot the absorption in relatively 
pure hydrogen for various plasma densities. At 
the lowest plasma densities Up =2.6 ma) at which 
the signal-to-noise ratio allowed a reproducible 
run, the a” line is at its cyclotron field and the 
H,* and H* lines are broad and already shifted 
towards lower B fields. As the plasma density 
increases, all three lines shift towards lower 
fields in qualitative agreement with Fig. 1. The 
shift from cyclotron resonance to hybrid reso- 
nance requires an increase in density of approxi- 
mately two orders of magnitude, in agreement 
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ABSORPTION IN RELATIVE UNITS 





























IHS 
fe) 7 14 21 28 35 42 ag 
MAGNETIC FIELD IN KILOGAUSS 





FIG. 2. Absorption and dispersion in a helium plas- 
ma with hydrogen impurity. Helium pressure= 100 
microns; hydrogen concentration unknown but less than 
2%. Quadrupole field, f=10.82 Mc/sec. Tube current 
=8ma. The arrows indicate the cyclotron fields of the 
various ions. 


with Fig. 1. A quantitative comparison with theory 
cannot yet be made. We do not know the relative 
concentrations of the various ions in the plasma 
which, moreover, are apt to vary as the plasma 
density is increased. There exists, however, one 
notable disagreement with Fig. 1. The H,* line 
definitely crosses the H,* cyclotron field and the 
H,* line crosses the H* field. The reason for 

this is not yet known, but it is thought to result 
from the presence of density gradients in the 
plasma. The behavior of linewidths with increas- 
ing density is of interest. The width of the H,* 


line increases, while that of the H,* line decreases. 


This is reasonable, however, for the following 
reasons. In the transition region between cyclo- 
tron and hybrid resonances, the frequency (and 
thus the resonant magnetic field) is density de- 
pendent. In a nonuniform plasma one then ex- 
pects the line to be additionally broadened. Now, 
in Fig. 3 the H,* line is in the low-density part 
of the transition region, while the H,* line is in 
the high-density part of it; it is this fact which, it 
is now thought, may account for the apparently 
anomalous behavior of the linewidths. 

Results similar to, but not identical with, those 
shown in Fig. 3 were obtained at other frequencies 
and pressures and for various other gases.° 

The author is grateful to H. W. Dail, N. R. 
Sheeley, and A. S. Wilczek whose aid made this 
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FIG. 3. Absorption in a hydrogen plasma for various 
plasma densities (tube currents), pressure=110 mi- 
crons. Dipole field, f=18.18 Mc/sec. The arrows 
indicate the cyclotron fields of the various ions. The 
origin of the structure in the absorption lines at high 
magnetic fields is not known. 


experiment possible. Many stimulating discus- 
sions with J. K. Galt are acknowledged. 





'E. Astrom, Arkiv Fysik 2, 443 (1950). 

2p, L. Auer, H. Hurwitz, Jr., and R. D. Miller, 
Phys. Fluids 1, 501 (1958). 

3K, Korper, Z. Naturforsch. al2, 815 (1957); a12, 
220 (1960). 

‘s. J. Buchsbaum, Phys. Fluids 3, 418 (1960). 

5s, J. Buchsbaum (to be published). 

°T, H. Stix, Phys. Rev. 106, 1146 (1957); T. H. 
Stix and R, W. Palladino, Phys. Fluids 1, 446 (1958). 
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CONTRIBUTION OF THE FERMI CONTACT TERM TO THE MAGNETIC FIELD 
AT THE NUCLEUS IN FERROMAGNETS* 


A. J. Freeman, 
Materials Research Laboratory, Ordnance Materials Research Office, Watertown, Massachusetts and Solid State 
and Molecular Theory Group, Massachusetts Institute of Technology, Cambridge, Massachusetts 


and 


R. E. Watson, 
AVCO, RAD, Wilmington, Massachusetts and Solid State and Molecular Theory Group, Massachusetts Institute of 
Technology, Cambridge, Massachusetts 
(Received November 10, 1960) 


Recent Méssbauer measurements’ have revealed constant (~ -3 atomic units) and this has led to 
a large but unexpectedly’ negative (-333 kgauss) the general assumption of spin polarization effects 
effective magnetic field, Hg, at the Fe*’ nucleus linear in total spin. 
in ferromagnetic iron. These results led to the We have calculated x for the free ions Mn*, 
recognition’ that the dominant source of Hz must Fe**, and Nit?; the resulting x’s are -3.34, -3.29, 
come from the polarization of the core s electrons and -3.94a.u., respectively. We also obtained a 
which then contribute a field, H,, through the x of -3.27 a.u. for Ni* in a crude (i.e., point 
Fermi contact interaction. We are here report- charge) crystalline field. The variation in y for 
ing results for a set of calculations* for H, which _ the Ni** calculations is due to its extreme sensi- 
are based on a model consistent with neutron dif- tivity to the behavior of the 3d spin distribution. 
fraction* and energy band results.° These results The calculated y’s are consistently more negative 
are compared with those just reported by Goodings than experimental values and this overestimation 


and Heine® for a less realistic model. Both cal- appears to be a general result of such calculations’ 
culations are shown to be deficient in that the ex- Marshall’s? estimates for the outer electrons in 
perimental value of H, is not realized, and hence metallic Co gave a contribution to the magnetic 
that the dominance of the core contribution has field of about +260 kgauss. Such estimates are 


not been proven. Some striking results are also necessarily crude, but probably 200 kgauss for 
reported for an investigation of H, as a function these terms in Fe is reasonable. (Goodings and 
of the position of the 3d density in the atom. Heine estimated this to be between +120 and +300 
These results help to clarify the theoretical pic- kgauss.) To give a net H, of -333 kgauss in Fe”, 
ture and should lead to a proper theoretical treat- the actual H, must be ~-530 kgauss (or there must 
ment of the problem. be another, as yet uninvestigated negative contri- 
Our investigations are based on “spin (exchange) bution). If we estimate by assuming, as for the 
polarized” Hartree-Fock® calculations, i.e., cal- divalent ions, that x is about -3 a.u., then H, for 


culations in which electrons in the same shell, the metal core (with 2.2 unpaired spins) is only 
but differing in spin direction (mg) are allowed -280 kgauss and this is clearly not enough. 
different radial distributions. Physically, this To account for this large difference between 
difference arises, in atoms having a net spin, core polarizations in metals and ions, several 
from the different exchange interactions experi- factors were considered. It is generally assumed 
enced by electrons of majority (say up, fT) and that in metallic iron there are between 7 and 8 
minority (1) spins. The resultant difference in electrons in the 3d band (as against 6 for the neu- 
radial distributions leads to a difference in den- tral atom or divalent ion). Furthermore, theore- 
sity for s electrons at the nucleus, | p,(0) - tical evidence® indicates that the 3d electrons in 
- |¥(0)1?, which contributes to Hg through the the metal are more expanded than in the free 
Fermi contact interaction. A convenient meas- (3d° 4s*) atom. We have therefore done a full 
ure of this interaction is given by Hartree-Fock spin (exchange) polarized calcula- 
tion for Fe in both 3d° 4s? and 3d® configurations. 
x= (40/3) erst! ¥4 |? - |p, (0) 17}, As shown in Fig. 1, even though the 3d® one-elec- 


tron 3d charge density is expanded relative to the 
where S designates the number of unpaired elec- 3d° 4s’, the resultant one-electron 3d spin density 
trons. For divalent transition element ions, ex- (the relevant factor determining H,)° is not signifi: 
perimental values of x are known’ to be roughly cantly changed. This spin density agrees with 
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, FIG. 1. One-electron 3d spin and charge densities for Fe 3d° 4s* and 3d°. 
d neutron measurements‘ and energy band calcu- their lower limit (+120 kgauss) to make the cal- 
. lations® which give a spin density which is actually culated H,. explain the experimental He. 
gi. slightly contracted relative to the free 3d° 4s? We do not agree. If their free-atom overesti- 
hust atom. For Ni*?, such a contraction in spin den- mate of H,, is corrected (from -355 to ~ -280 
ad sity was shown to reduce H,..° kgauss), then their expanded 3d value of H, would 
‘ For the 3d° state, the calculated H; (-350 kgauss) be changed from -420 kgauss approximately back 
” is 30 kgauss lower than the 3d°4s’ value (both are _—to their quoted free-atom value. Furthermore, 
| overestimated by about 40 kgauss). This improve- their “artificial 3d expansion” gives a spin den- 
ment is too small to overcome the positive contri- sity in direct disagreement with the neutron dif- 
butions and so cannot account for the experimental fraction and energy band results cited above. 
result. An investigation of H, as a function of the posi- 
oa Goodings and Heine® did a 3d°4s* free-atom spin _ tion of the 3d density in the atom helps to clarify 
polarized calculation for Fe and obtained a x of the theoretical picture. For a series of fixed 3d 
- -3.8 a.u. (-355 kgauss) from the core. The dif- functions (which for computational simplicity we 
- ference between this value and our own arises took as a single Slater atomic orbital), we deter- 
7 from the computational sensitivity of the calcu- mined H, via spin-polarized Hartree-Fock calcu- 
lations (thus the best one can hope to do is to lations. The results are shown in Fig. 2 asa 
provide a consistent set of computations). Their function of the maxima of the 3d orbitals. 
o x (and ours) shows the overestimation discussed We see that a large, but unrealistic, expansion 
P above. To estimate the effect in a metal, they or contraction of the 3d shell spin density will 
» repeated the calculation with fixed “3d functions yield an H,, which is more than sufficient to “ex- 
Me artificially expanded by about 5% at the maxima plain” the metal results. Simple scaling does not 
ity and 10% over the tail regions.” This lowered appear to be the answer to the problem.® On the 
al H, (by 65 kgauss) to -420 kgauss. This, they felt, other hand, Fig. 2 suggests that any (other) ad 
sufficiently balanced the positive contribution at hoc small change of shape in the 3d functions 
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FIG. 2. Core effective field, H,, as a function of 
the maximum of the 3d orbital for Fe 3d°. 


which removes spin density from the region rough- 
ly between 0.5 and 0.75 a.u. will yield a substantial 
lowering of H,. This can be devised while main- 
taining reasonable compatibility with the neutron 
measurements. 

What is needed is a theoretical treatment which 
gives a more realistic description of the atom’s 


environment, the importance of which is indicated 
by the experimental results’® for Fe*”’ in Fe, Co, 
and Ni. 

We are indebted to W. Marshall, R. K. Nesbet, 
and R. J. Harrison for valuable discussions and 
assistance. 





*Some of the material in this Letter was discussed 
at the Sixth Conference on Magnetism and Magnetic 
Materials, November 14-17, New York City, Suppl. J. 
Appl. Phys. (to be published). 

+The work of this author has been supported in part 
by the Ordnance Materials Research Office. 
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"See a recent review by G. K, Wertheim, Suppl. J. 
Appl. Phys. (to be published). 





OPTICAL ACTIVITY IN TELLURIUM 


K. C. Nomura 
Honeywell Research Center, Hopkins, Minnesota 
(Received October 21, 1960) 


The crystal structure’ of the elemental semi- 
conductor, tellurium, has a symmetry which 
places it in an enantiomorphous class. Since this 
kind of symmetry constitutes only a necessary 
condition for the existence of optical activity, 
this effect was examined in single crystals and 
tellurium was found to have a large rotatory 
power. Thus, this is the first of the space-group- 
dependent properties to be measured in this semi- 
conductor. 

Both left- and right-handed single crystals were 
grown by the Czochralski method? which corres- 
pond to the space groups D,° and D,*, respectively. 
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Such crystals are levo- and dextrorotatory. At 
the initial stages of this work, all of our crystals 
were found to be dextrorotatory. This, of course, 
was to be expected since all of our crystals were 
“descended” from the same parent seed of a few 
years ago. The levorotatory crystals were grown 
from seeds which were randomly grown as nee- 
dles in a hydrogen atmosphere. We would reason- 
ably expect that such needles would have equal 
likelihood of being left- or right-handed. Quite 
fortuitously, a levorotatory crystal was grown 
on the second attempt. 

The samples were cut into thicknesses varying 
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from 0.212 cm to 1.12 cm with plane parallel 
faces on the (0001) plane. The faces were lapped, 
polished in chromic acid, and oriented on (0001) 
with an accuracy of better than 0.5 degree. The 
latter is particularly important to prevent spuri- 
ous interference effects.* The measurements 
were made using a Perkin-Elmer, Model 112, 
single-beam double-pass infrared spectrometer 
with NaCl optics. Stacked plates of AgCl were 
used for polarizers, and advantage was taken of 
the inherent 30% instrument polarization by plac- 
ing the second polarizer in the same polarization 
orientation. A well-collimated plane polarized 
beam was directed normal to the first face of the 
crystal, and the sample was rotated to insure that 
misorientations and misalignment effects were 
not present. This general procedure was found 

to give results which were in agreement with the 
published values* for natural single crystals of 

a quartz. 

It is conventional in optical activity experiments 
to express the specific rotatory power (@/d) 
(quoted in units of degrees of rotation per mm of 
sample thickness) as a function of wavelength or 
wave number. The results for both levo- and 
dextrorotatory crystals are shown in Fig. 1, in 
which the abscissa is expressed in terms of the 
wave number ? (in cm™*) so that deviation from 
normal rotatory dispersion can be seen. 

These results may now be analyzed in terms of 
Fresnel’s expression for the specific rotatory 
power: 


6/d =18(n Ro” yp” (degrees /mm) 


where 7 is incm™. The quantity, np-nz, is the 
difference in refractive indices for the ordinary 
and extraordinary rays in the direction of the optic 
axis. At )=0.20x10* cm™ (5 microns), np -n7 
=+1.58x107%, values some fifty times greater 

than for a quartz at A =0.589 micron. 

From the curves of Fig. 1, it is seen that de- 
viations from normal dispersion occur for 7 >0.20 
x10*cm™. This is consistent with refractive 
index results reported by Hartig and Loferski® 
and by Caldwell and Fan.® 

It is a pleasure to thank L. P. Beardsley for 
his capable assistance, A. Ryberg who grew the 
crystals, and J. S. Blakemore for many useful 
discussions. 
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FIG. 1. The rotatory powers of dexro- and levo- 
rotatory tellurium as a function of wave number. 
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ROOM-TEMPERATURE DISLOCATION DECORATION INSIDE LARGE CRYSTALS* 


Charles B. Childs 
National Aeronautics and Space Administration, Goddard Space Flight Center, Greenbelt, Maryland 


and 


Lawrence Slifkin 
University of North Carolina, Chapel Hill, North Carolina 
(Received November 8, 1960) 


The usual techniques for decoration of disloca- 
tions in ionic crystals! generally either require 
a high-temperature annealing treatment or else 
display only those dislocations near the crystal 
surface. For the study of some aspects of the 
plasticity of ionic crystals, it would be desirable 
to have a means of observing dislocation configu- 
rations throughout large regions of the interior 
of the specimen without the necessity of an anneal- 
ing procedure which also changes the dislocation 
arrays. This note describes such a technique 
which displays with great contrast the distribution 
in the interior of large single crystals of silver 
chloride. 

The work of Castle? has demonstrated that if 
internal print-out silver is deposited within the 
crystal by means of sweeping in surface-produced 
photoelectrons (the Haynes-Shockley method), one 
can see something of the dislocation structure of 
the crystal. Since direct observation was pre- 
sumably difficult because of the excessive amount 
of light scattered by random, submicroscopic 
silver specks, Castle sectioned his specimen and 
looked at the microscopic silver precipitates re- 
vealed on each freshly cut surface. This method 
thus did not yield the large-region views produced 





FIG. 1. A square net of dislocations. 
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by the precipitation processes of Amelinckx' and 
of Mitchell and co-workers.' 

The essential contribution of the present note is 
the discovery that, after production of only moder- 
ate print-out densities by means of the Haynes- 
Shockley method, a room temperature aging of 
the crystal results in extensive migration of silver 
from random submicroscopic specks into the dis- 
locations. After a period of several days the dis- 
locations have become so sharply decorated re- 
lative to the diffuse background that they may be 
plainly observed even with ordinary bright field 
microscopy and a simple 20x objective, the work- 
ing distance of which is great enough to allow of 
continuous examination deep into the crystal. In 
addition to random dislocations, structures such 
as pure tilt boundaries, square and hexagonal nets, 
and helices have been observed. The accompany- 
ing photographs (Figs. 1 and 2) illustrate the reso- 
lution and clear contrast achieved. The five speci- 
mens employed were Harshaw crystals which had 
previously been annealed in air on quartz powder. 

It appears, then, that whereas the electron 
traps of greatest total cross section inside the 
crystal are not dislocations (perhaps they are 
impurity atoms), the most stable sites for the 
final precipitated silver are largely at disloca- 
tions. The migration of silver is most probably 
as ions plus electrons. Suptitz®’ has shown that 
the rate of this migration is greatly enhanced by 





FIG. 2. A hexagonal net of dislocations. 
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irradiation of the crystal with visible light, pre- 
sumably a photoionization effect. This suggests 
the possibility of accelerating the decoration pro- 
cess by long-wavelength irradiation. It would 

also appear that any effect of internal print-out 

on the plastic yield stress of silver chloride crys- 
tals would initially be relatively small, but would 
increase greatly as the silver migrated to the dis- 
locations. Miller,‘ in this laboratory, has indeed 
found that the initial hardening is only slight; ex- 
periments are now under way to investigate the 
effects of aging. 

Among the experiments suggested by this deco- 
ration technique are studies of possible internal 
dislocation sources and stress-induced pileups 
in silver chloride. Since any slow room tempera- 


ture recovery is expected to be arrested at an 
early stage in the decoration process, it should 
be possible to study the dislocation distributions 
resulting from plastic deformation without the 
necessity of an intervening opportunity for their 
rearrangement. . 





*This research was performed at the University of 
North Carolina. The support of the University Research 
Council and the Air Force Office of Scientific Research 
is gratefully acknowledged. 

‘See, for example, J. T. Bartlett and J. W. Mitchell, 
Phil. Mag. 5, 445 (1960); and S. Amelinckx, Phil. Mag. 
3, 307 (1958). 

2J. W. Castle, J. Appl. Phys. 28, 743 (1957). 

3p, Suptitz, Z. wiss. Photographie 53, 201 (1959). 

‘M. Miller (to be published). 





ELECTRONIC CONDUCTION AND EXCHANGE INTERACTION IN A NEW CLASS 
OF CONDUCTIVE ORGANIC SOLIDS 


R. G. Kepler, P. E. Bierstedt, and R. E. Merrifield 
Central Research Department, Experimental Station, 
E. I. du Pont de Nemours and Company, Wilmington, Delaware 
(Received November 7, 1960) 


A group of highly conductive organic solids 
which offer a unique opportunity for experimental 
study of electronic conduction and interaction in 
solids has recently been discovered in this labo- 
ratory. These materials are salts of the radical- 
anion formed by addition of an electron to tetra- 
cyanoquinodimethan (TCNQ): 


(NC),C = <>» =C(CN),. 


TCNQ 


In combination with a variety of cations, the 
[TCNQ:]  radical-anion forms two series of salts 
having the compositions 


M*[TCNQ-] and M*[TCNQ-] [TCNQ]°, 


respectively.’ In both types of salt the band 
formed from the lowest normally unfilled mole- 
cular orbital of the TCNQ is partially filled, 
being half filled for the first series of salts and 
one-quarter filled for the second. As a result of 
the relatively weak intermolecular interactions 
characteristic of molecular crystals, this band 
will be quite narrow. In addition, both the band- 
width and the exchange coupling of the odd elec- 
trons* appear to be strongly influenced by the na- 


ture of the cation M*, and by whether or not 
[TCNQ]° is present in the crystal. 

In the quinolinium salt (see Table I), the elec- 
trons of the [TCNQ-] radicals appear to form a 
conventional degenerate system similar to a 
metal. This is evidenced by the relatively large 
electrical conductivity (the largest yet reported 
for an organic solid), the virtual absence of 
activation energy for conductivity, and the tem- 
perature-independent paramagnetism. The mag- 
nitude of the paramagnetic component of the 
susceptibility corresponds to a density of states 
at the Fermi surface of 7.5 states/ev molecule, 
which implies a quite narrow energy band. 

In the triethylammonium salt the paramagnetic 
component of the magnetic susceptibility is given 
by 


Xp= (2g*6°N/kT)[3+exp(J/kT)}", (1) 


with J=0.041 ev. This corresponds to an assem- 
bly of N quasi-molecules, each having a triplet 
state lying at an energy J above a singlet ground 
state. A single sharp absorption line is observed 
in the electron spin resonance spectrum at 
g=2.002.2 The temperature dependence of the 
intensity of this line also follows Eq. (1). The 
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Table I. Electronic properties of representative radical-anion salts. 
Electrical conductivity Magnetic susceptibility 
Value Activation Value Temperature 
Material at 300°K energy at 300°K dependence 
(ohm~ cm-! (ev) (emu-mole™') 
O\ 
(TCNQ:] [TCNQ]® 1004 <0.01 +2.2x10+4 Decreases gradually to 
\ \y +1.0x10~ at 77°K 
i 
+ o 0 a ~? 
[(C,Hs);NH]  [TCNQ-] [TCNQ] 4.0 0.14 +6.4x10 Eq. (1) with J=0.041 ev 
- c 
K* [TCNQ:] 1.9x10+ 0.36 -1.3x104 Temperature independ- 


ent from 77° to 450°K4 





“In direction of highest conductivity (see text). 
bparamagnetic contribution. 
°Measured on compressed powder sample. 


dj~0.15 ev from temperature dependence of spin resonance. 


electrical conductivity of this salt is highly ani- 
sotropic, having room-temperature values of 
4.0, 0.05, and 0.001 ohm~'cm™ in the three 
principal crystal directions and varies exponen- 
tially with temperature with an activation energy 
of 0.14 ev, which is isotropic within experimen- 
tal error. Preliminary crystallographic results 
show that this salt belongs to the monoclinic C2 
space group with four formula units per unit cell. 
The TCNQ units appear to be arranged in infinite 
face-to-face stacks with the direction of highest 
conductivity lying along these stacks and approxi- 
mately normal to the planes of the TCNQ mole- 
cules. The thermoelectric power is -100 yv/°C 
in the direction of highest conductivity, which 
indicates that electrons are the majority car- 
riers in this material. We have been unable to 
detect a Hall voltage in this salt and estimate 
that such a voltage would have been observed if 
the electron mobility had been greater than 0.04 
cm?/v-sec.* A lower limit for the electron mo- 
bility of 0.02 cm?/v-sec can be calculated from 
the observed conductivity and the assumption that 
all of the electrons of the radical-anions partici- 
pate in carrying current. 
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The final example in Table I, K* [TCNQ-] , pro- 
vides a case in which there is a very low electri- 
cal conductivity and a relatively high activation 
energy for conductivity. The magnetic suscepti- 
bility is negative and temperature dependent al- 
though there are indications from spin resonance 
measurements of a triplet level lying ~0.15 ev 
above the ground state. 

We are currently studying a wide range of 
solids containing [TCNQ:] radical-anions and will 
report on these in detail elsewhere. 

We wish to acknowledge the assistance of Dr. 
Paul Arthur, Jr. in providing us with preliminary 
results from his crystallographic studies on some 
of these materials. 





'D. S. Acker et al., J. Am. Chem. Soc. (to be pub- 
lished). 

*D. B. Chesnut, H. Foster, and W. D. Phillips, J. 
Chem. Phys. (to be published). 

3It is possible that charge carrier motion in these 
materials takes place by an activated hopping process, 
in which case it is questionable whether a Hall effect 
should be observed; see discussion in J. Phys. Chem. 
Solids 8, 50 (1959). 
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DIRECT OBSERVATION OF EXCITON MOTION IN Cast 


D. G. Thomas 
Bell Telephone Laboratories, Murray Hill, New Jersey 


and 


: J. J. Hopfield — 
Bell Telephone Laboratories, Murray Hill, New Jersey and Ecole Normale Supérieure, Paris, France 
(Received November 7, 1960) 


Several authors have recently described or pre- 
dicted a number of rather unusual optical effects 
connected with the exciton spectra of crystals. 
Gross and Kaplyanskii’ have shown that the weak 
forbidden transition to the n=1 state of an exciton 
in cubic Cu,O has absorption and polarization 
properties which depend upon the direction of the 
exciting light in the crystal. Pekar and Tsekvava? 
have predicted that cubic crystals can show bire- 
fringence near exciton absorption peaks if the 
exciton has an anisotropic mass. In CdS oriented 
with the hexagonal c axis perpendicular to a mag- 
netic field H, and with the light traveling at right 
angles to both, Hopfield and Thomas’ found that 
the exciton spectra depended markedly on the 
sign of H. These effects are, however, not pe- 
culiar to excitons, since quadrupole transitions 
of atoms in a crystal field can exhibit analogous 
behavior. Some effects peculiar to excitons are 
the result of excitons having an energy versus 
crystal momentum curve which is not flat. Thus 
Pekar* has shown theoretically that for wave- 
lengths close to a suitable exciton transition, 
light of the same frequency and polarization can 
be propagated in the same direction with different 
refractive indices; Brodin and Pekar® have offered 
evidence of this effect in anthracene crystals at 
20°K. In this Letter we report on a new experi- 
ment which gives a direct observation of the curva- 
ture of the E(k) relation for excitons. 

In reference 3 it was pointed out that an intrinsic 
direct exciton created from a photon must have 
the same crystal momentum vector as the photon. 
As a result of its momentum the exciton has a 
velocity Ak/M, where k is the wave vector of the 
light in the medium and M is the exciton mass. 

In a magnetic field this velocity gives rise to an 
additional term in the exciton Hamiltonian equi- 
valent to a quasi-electric field Eg = hk xH/Mc. 

Stark effects on excitons are in general even func- 
tions of the electric field. This will be true for 

an applied electric field for H=0. However, in 

the presence of a magnetic field the effective elec - 
tric field is the sum of the applied and quasi-fields. 
When the applied electric field, the magnetic field, 


and k are all mutually perpendicular the effect of 
the quasi-field should be to shift the plane of sym- 
metry of the Stark effect by an amount Eg. From 
Eq and k the exciton mass can be directly calcu- 
lated. 

The Stark effect has been observed on the n =2 
exciton states derived from the top valence band 
of CdS at 1.6°K in light polarized paralled to c. 
The Stark effect of the n=2 states is used since 
it is much larger than the effect expected for the 
n=1 states. CdS is piezoelectric and so an elec- 
tric field could affect the spectra by changing the 
band gap; however, for all fields used no change 
in the positions of the n=1 lines could be detected, 
indicating that the piezoelectric effects are neg- 
ligible. The crystals were 5-10 microns thick 
and several mm long; the electric field was ap- 
plied by placing the crystals between the plates 
of a condenser to which a symmetrical square 
wave was applied, thus avoiding the necessity of 
making contacts. The crystals were illuminated 
with nearly monochromatic light from a Perkin 
Elmer monochromator with a mean wavelength 
equal to that of the exciton transition. This en- 
sured that the crystal was exposed to the mini- 
mum light intensity so that the photoconductivity 
was as small as possible and the effective dielec- 
tric relaxation time as long as possible. The 
transmitted light was photographed with a high- 
dispersion Bausch and Lomb grating spectrograph; 
in front of the entrance slit there was a revolving 
sectored disk which served the dual purpose of 
switching the square wave by means of a separate 
light source and photomultiplier, and of exposing 
the photographic plate during only one phase of 
the square wave. Hence a spectrum could be ob- 
served for the electric field in either of the two 
possible directions. 

Figure 1 shows the effect of the electric field 
on the separation of two pairs of lines of the Zee- 
man pattern of the n =2 exciton state in CdS with 
c1H and H = 31000 gauss [the symbols are identi- 
fied in reference 3]. The two curves are for the 
two directions of H which, as explained above, 
give different spectra. This orientation of H with 
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FIG. 1. The separation between two pairs of lines in 
the Zeeman spectra of the n=2 exciton state in CdS at 
1.6°K, as a function of the applied electric field. The 
pair of lines P,~, P,~ are two lines in the absorption 
spectrum with H in one direction, and the pair of lines 
P,*, Py* are two analogous lines for H in the other 
direction. The light is polarized parallel to the hexa- 
gonal c axis. 


respect to c was chosen as it gave the most 
clearly defined lines. The duration of the square- 
wave cycle was 400 usec in these experiments. 

In order to show that the frequency is high enough 
to prevent appreciable screening of the electric 
field by polarization of the crystal, a square- 
wave generator was used without chopping the 
light. It was found that for H =0 there was no 
change in the effect as the duration was changed 
from 1.6 milliseconds to 20 psec. It is seen that 
the two curves are approximately symmetrical 
with their planes of symmetry on either side of 
the plane of zero applied field. In the absence of 
a magnetic field no such asymmetry was found. 
The mean displacement of the planes of symmetry 
from the center point is about 95+ 10 v/cm. The 
directions of the displacements were found to be 
in accord with the theory described above. At 
high applied electric fields the lines become 
broadened and there is evidence that the field in 
the crystal is reduced. If Eg equals 95 v/em 
then M is found to be 1.1mg. (The refractive 
index for E ||c at the frequency of the n=2 tran- 
sition is estimated from the reflectance spectra 
to be 2.2; this value was used to calculate the 
photon momentum inside the crystal.) A recent 
analysis® of the exciton Zeeman spectra of CdS 
has yielded an exciton mass in a direction per- 
pendicular to c equal to0.9m,, where mg is the 
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free electron mass. The observation of the effect 
and its agreement in sign and approximate mag- 
nitude with the predicted values demonstrate 

the exciton momentum and show that this exciton 
“mass spectroscopy” experiment can be used to 
determine exciton masses. 

In some crystals the Stark effect was found to 
be considerably larger than shown in Fig. 1 and 
the apparent value of Eg considerably less, in 
one case only 30 v/cm, although with the correct 
sign. These effects can be understood by sup- 
posing that at 1.6°K and 31000 gauss in some 
crystals tangent (Hall angle) >1, i.e., that w-7>1, 
where wc is the cyclotron resonance frequency 
of the highest mobility carriers (presumably elec- 
trons) and 7 is the relaxation time. Under these 
circumstances a transverse Hall field exceeding 
the applied field can form across the thin dimen- 
sion of the crystal in the direction of the light 
propagation. This field can be formed before the 
longitudinal field is screened out because screen- 
ing of the longitudinal field by surface charge is 
much less efficient than that of the transverse 
field. Presumably at sufficiently high frequencies 
the Hall field would not be established, but such 
frequencies were not available to us. In one crys- 
tal, analysis of the results showed that w-7~ 1.7. 
This indicates that cyclotron resonance may be 
observed in CdS but only in selected crystals at 
fields in excess of about 20000 gauss. From the 
magnitude of the Stark shift shown in Fig. 1, it 
is estimated that w-7~0.5 for this crystal and 
that the Hall field shifts the position of the sym- 
metry planes by about 20%. The calculated ex- 
citon mass is therefore expected to be somewhat 
too high. 

The method of determining the exciton mass 
described here could be of use in determining 
the individual masses of holes and electrons since 
the reduced exciton mass comes from the exciton 
ionization energy. Interest also derives, how- 
ever, from the fact that due to polaron and other 
effects the exciton mass may not be simply the 
sum of electron and hole masses. Provided the 
exciton lines are sharp, which was the case for 
the crystal used in Fig. 1, the method is easiest 
to apply when WeT<1. Finally, the inclusion of 
the Ak xH/Mc term in the Hamiltonian is of im- 
portance in analyzing the Zeeman effect of ex- 
citons, particularly for the higher quantum num- 
ber exciton states. 

Thanks are due to Dr. D. C. Reynolds of Wright 
Patterson Air Force Base for some of the crys- 
tals used in this work, and to E. A. Sadowski for 
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PROTON MOTIONS IN AMMONIUM HALIDES BY SLOW 
NEUTRON CROSS-SECTION MEASUREMENTS* 


J. J. Rush, T. lL. Taylor, and W. W. Havens, Jr. 
Departments of Chemistry and Physics, Columbia University, New York, New York 
and Brookhaven National Laboratory, Upton, New York 
(Received November 8, 1960) 


The motions of ammonium ions and the order- 
disorder transitions in their salts have been in- 
vestigated extensively by a variety of techniques. 
These include specific heat measurements,'~* 
crystal structure studies by neutron diffraction,° 
infrared spectroscopy,®~® and nuclear magnetic 
resonance.*-" The results of the investigations 
reported here show that cross-section measure- 
ments for neutrons of subthermal energy also 
provide a valuable criterion for determining the 
type of proton motion and the strength of chemical 
binding in the ammonium halides or other com- 
pounds of hydrogen. 

Using a crystal spectrometer™,** at the Brook- 
haven reactor, total neutron cross-section meas- 
urements were made at room temperature for 
NH, gas at 115 psia and for the following ammoni- 
um halide salts: NH,F, NH,Cl, NH,Br, and NH,I. 
In Fig. 1 the cross sections per hydrogen atom, 
corrected for neutron absorption,* are plotted 
as a function of neutron wavelength. Since the 
scattering cross section for hydrogen is large 
and primarily incoherent, Bragg interference 
effects are relatively small and do not contribute 
significantly to the total cross section. A series 
of closely spaced cross-section measurements 
in the neighborhood of the most intense Bragg 
peak for NH,Br indicated that the Bragg effects 
were small and could be neglected. 

At short wavelengths (),,~ 0) the scattering 
cross sections approach the values for the free 
atoms. At long wavelengths the predominant in- 
elastic scattering process is one in which the 


neutron gains energy from the crystal or molecule. 


The scattering cross section then becomes pro- 
portional to the wavelength, A,, and may be re- 


presented by 0, =a+bd,. The slopes, 0, of the 
cross-section curves in Fig. 1 were determined 
in the long-wavelength range by a least-squares 
fit to the data. The dashed line is the theoretical 
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FIG. 1. Neutron scattering cross section per hydro- 
gen atom for ammonia and the ammonium halides. 
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cross section per proton in NH, calculated using 
the Krieger-Nelkin method.** 

The slopes from Fig. 1 together with other pro- 
perties of the ammonium halides are summarized 
in Table I. Previous work referred to above has 
established that the ammonium ions in phases II 
and III may be considered torsional oscillators 
with the frequencies listed and with barriers to 
rotation® well above kT,. The transition from 
phase III to phase II at T) is a cooperative order- 
disorder process involving reorientation of the 
ammonium ions. At room temperature NH,]I is 
in phase I where the ammonium ions are thought 
to be freely rotating around the N-H --I™ axes 
with a potential barrier to rotation of approxi- 
mately 35 cm™?.” 

The slopes and the absolute magnitudes of the 
scattering cross sections show a definite corre- 
lation with the proton dynamics and hydrogen 
bonding in the ammonium halide salts. As shown 
in Table I, the slopes decrease with increasing 
torsional frequencies, although a more detailed 
comparison would involve consideration of the 
entire fundamental and lattice vibrational spec- 
trum. The order of the torsional frequencies 
and the slopes is that expected from the abilities 
of the halide ions to form hydrogen bonds, i.e., 
F'>Cl°>Br’ >I. Thus, the strong hydrogen 


bonds in ammonium fluoride result in a relatively 
high torsional frequency with fewer states ex- 
cited and a low slope (2.85 barns/A). The weaker 
bonds in the other ammonium halides and a struc- 
ture which allows reorientation of the ammonium 
ions result in higher slopes. 

The slope of the scattering cross-section curve 
for NH,I at room temperature approaches that 
for NH, gas. This is consistent with the model 
in which the ammonium ion in phase I of NH,]I is 
relatively free to rotate about one axis, as estab- 
lished from other evidence. Preliminary meas- 
urements of the slope of the cross-section curve 
for NH,I at 0°C and -20°C result in values of 11 
and 6 barns/A, respectively. Thus, there is a 
marked change in slope in passing from phase I 
in which the ammonium ion is freely rotating to 
phase II in which it is a torsional oscillator. Re- 
latively little change was observed in passing 
through the order-disorder transition from phase 
III to phase II for NH,Cl. It is interesting to note 
that the slopes for NH,Br and NH,Cl are meas- 
urably different even though the difference in 
their torsional frequencies is small. The limit- 
ing slope of the scattering cross-section curve 
appears to be correlated with the number of ex- 
cited states and therefore the freedom of motion 
of the hydrogen atoms. 


Table I. Limiting slopes of the scattering cross-section curves for the ammonium halides compared with other 
properties. The errors listed include those due to counting and to uncertainty in the absorption cross sections. * 








NH,F NH,Cl NH,Br NH,I 

Slope, barns/A, 23°C 2.85 +0.21 4.83 +0.18 5.71 +0.09 11.21 +0.18 
Wavelength range, A 8-11 6-11 6-11 7-11 
Cross section per proton 

at 10 A, barns 105 125 128 162 
Motion of NH,* at 23°C ose. ose. osc. rot. 
Transition temperatures, °C 

Phase I (NaCl) +> Phase II (CsCl) eee 184.3 137.8 -17.6 

Phase II (CsCl) -—~ Phase III» tee -30.5 -38.1 -41.6 
Torsional frequency, cm~! © 

Phase III 5235 391 319 279 

Phase II eee 359 311 eee 
Potential barrier to rotation® 

kcal/mole (Phase II) tee 5.3 4.2 2.9 








8See reference 14. 


phase III for NH,Cl has the cubic CsCl structure, while that for NH,Br and NH,I is slightly distorted to a tetra- 
gonal structure. NH,F has a hexagonal structure in which each NH,* ion is surrounded tetrahedrally by four F~ 
ions. The only observed thermal transition is a small one about -30°C (references 1 and 7). 


“See references 6 and 7. 
dsee reference 9. 
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The data presented here show that total cross- 
section measurements can give valuable evidence 
concerning the dynamics and chemical binding 
of hydrogen in molecules and crystals. Meas- 
urements are now being extended to other am- 
monium salts and hydrogen compounds. 
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ROTATING SELF-CONSISTENT FIELDS AND ROTATIONAL STATES OF NUCLEI 


D. J. Thouless* and J. G. Valatin* 
Institute for Theoretical Physics, University of Copenhagen, Copenhagen, Denmark 
(Received October 7, 1960) 


Moments of inertia of rotating nuclei have been 
calculated by means of the cranking model’~* in 
which the particle motion is investigated in a po- 
tential rotating with a given angular velocity. 

The aim of the present work is to show that the 
self-consistent equations of the time-dependent 
Hartree-Fock method have such rotating solu- 
tions. The rotating potential is produced by the 
particles themselves and the particles move in 
this potential. Similar considerations hold for 
the generalized Hartree-Fock equations which 
take into account the pairing correlations. 

The time-dependent Hartree- Fock equations 
can be written as* 


iap/at=[v, pl, (1a) 
where p is the density matrix which satisfies the 
condition p?=p, and pv is the self-consistent en- 
ergy. The matrix elements of v are given by 
1 
Os aw?! 
v(x ,x on V75(x - x’) 
+ fdx,dx,'[(xx,|V |x’x,’) - (xx, |V 1x,'x’)]p(x,’,x,). 


(1b) 


For a system with translational invariance, the 


solutions of the time-independent equation 
[v,p] =0 (1c) 


are degenerate, because the center of the self- 
consistent potential may be moved to any posi- 
tion. This degeneracy permits solutions of the 
time-dependent equation (la) to be constructed, 
which represent a translational motion of the 
system with constant velocity.° The simplest 
way to obtain these solutions is to transform the 
equations to a reference system moving with a 
constant velocity v. Since there are no inertial 
forces acting in this system, the stationary so- 
lutions in the moving system can be obtained ex- 
plicitly in terms of the solutions with v =0. 

For a nucleus with a nonspherical equilibrium 
shape, the solutions of (1c) are degenerate be- 
cause of an arbitrary orientation of the axes of 
the system. The time-dependent equation (1a) 
has accordingly solutions corresponding to a ro- 
tation of these axes. They can be investigated by 
transforming the equation to a rotating reference 
system. Time-independent solutions in the ro- 
tating frame represent rotating solutions in the 
rest system, except in the case of rotation about 
a symmetry axis. Transforming (la) to a refer- 
ence system rotating with uniform angular velo- 
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city 2 about the y axis, one obtains the equation 
idp/ eat = [b+al,,A], (2a) 


where Ll, is the angular momentum operator of a 
single particle and the notation (, ) indicates that 
these quantities are expressed in terms of var- 
iables of the rotating system. Solutions of the 
time-independent equation 


[b+ Ql, ,p]=0 (2b) 


will reflect the effect of the centrifugal and 

Coriolis forces represented by the term Qly. 
For small values of 2, one can expand the so- 

lutions of (2b) in powers of 2, by writing 


p=pM+pM+..., D=p%+ M+... (3a) 


into the equation. One has [p,)]=0, and the 
terms linear in 2 give 


[5 p™]+[9,5]=-[01,,6] (8) 


The expression for »™ is obtained from the sec- 
ond term of (1b) with p replaced by p™. Assum- 
ing p® and »™ as known, (3b) is a linear equa- 
tion to determine p™. The moment of inertia 
obtained from the expression 


(ly) ™=trace Lp” (3c) 


differs from that of the cranking-model formula 
with independent-particle energies given by the 
eigenvalues of > through the correction term in 
p™, This represents the change in the rotating 
self-consistent potential due to the inertial forces 
acting on the particles. 

In the representation in which jp, > are di- 
agonal, (3b) is identical with equations obtained 
previously® by minimizing the ground-state en- 
ergy of the Hartree-Fock method with a pre- 
scribed expectation value of a total angular mo- 
mentum component. Equation (2b) would indeed 
result by minimizing the expectation value of the 
Hamiltonian in a rotating reference system, and 
this has the same expression as the Hamiltonian 
of the rest system plus a Lagrangian multiplier 
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times an angular momentum operator. This 
variational problem must certainly have solu- 
tions, so steadily rotating solutions of Eq. (1a) 


must exist. 
In the generalized Hartree-Fock method’ the 


ground-state energy depends on the density 
matrix p(x,x’), and on the expectation value 
x(x,x’) of a pair of field operators. These two 
quantities define the components of a projection 
operator K which acts in a single-particle space 
of twice as many dimensions, formed by the 
single-particle states and their adjoints. In addi- 
tion to (1b), there is a self-consistent potential 
determined by x, closely related to the energy 
gap of the elementary excitations. The two quan- 
tities define components of a matrix 9, and the 
equations of the generalized Hartree- Fock me- 
thod can be written in the form® 


ix / at = fn,x]. (4) 


The rotating solutions of this equation result then 
in complete analogy to those of (1a). 





*On leave from the Department of Mathematical 
Physics, University of Birmingham, Birmingham, 
England. 

'D. R. Inglis, Phys. Rev. 96, 1059 (1954); 97, 701 
(1955); A. Bohr and B. Mottelson, Kgl. Danske Viden- 
skab. Selskab, Mat.fys. Medd. 30, No. 1 (1955); J. G. 
Valatin, Proc. Roy. Soc. (London) A238, 132 (1956). 

24. Bohr and B. Mottelson, Norske Videnskab. Sel- 
skab, Forhandl. 31, No. 12 (1958) (to be consulted 
also for further references). 

3s. T. Belyaev, Kgl. Danske Videnskab. Selskab, 
Mat. fys. Medd. 31, No. 11 (1959). 

‘Pp. A. M. Dirac, Proc. Cambridge Phil. Soc. 26, 
376 (1930). 

5Time-dependent Slater determinant solutions of the 
Hartree-Fock method describing the translational mo- 
tion of nuclei have been considered previously by 
C. van Winter (private communication). 

8p, J. Thouless, Nuclear Phys. (to be published). 

"J. Bardeen, L. N. Cooper, and J. R. Schrieffer, 
Phys. Rev. 108, 1175 (1957); J. G. Valatin, Nuovo 
cimento 7, 843 (1958); N. N. Bogolyubov, Nuovo cimen- 
to 7, 794 (1958); Soviet Phys. — Uspekhi 67(2), 236 (1959). 
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EFFECT OF NUCLEAR ROTATION ON THE PAIRING CORRELATION 


Ben R. Mottelson and J. G. Valatin* 
NORDITA and The Institute for Theoretical Physics, University of Copenhagen, Copenhagen, Denmark 
(Received November 8, 1960) 


The empirically observed energy gap, the re- 
duction in the effective moment of inertia below 
that of a rigid rotator, and the transition to 
spherical shape with the approach to closed 
shells all testify to the important role of the 
pairing correlation in the low excited states of 
deformed nuclei. The total binding energy asso- 
ciated with the pairing correlations is, however, 
quite small,’ and thus a rather slight disturbance 
of the system may be sufficient to destroy the 
pairing and cause a transition to a less corre- 
lated particle motion. 

An analogy with the correlations in supercon- 
ductors has proved to be rather fruitful in inter- 
preting the nuclear features mentioned above.’ 
Pursuing this analogy, the close formal corre- 
spondence of the equations of motion in a constant 
magnetic field and in a rotating reference sys- 
tem suggests that the critical magnetic field 
phenomena of superconductors should also have 
their counterpart in the rotational spectra of 
nuclei. Indeed, the pairing force couples two 
particles in time-reversed single-particle states, 
and a rotation has an opposite effect on the par- 
ticles forming the pair. The Coriolis forces act 
in opposite directions and tend to decouple the 
pairing correlations. This is in analogy with the 
effect of the magnetic forces in a metal, which 
act with a different sign on electrons moving in 
opposite directions and destroy the pairing when 
the field reaches a critical value. The observa- 
tion or nonobservation of this effect in the rota- 
tional spectra of nuclei could serve as a test of 
the description of the pairing correlations in nu- 
clei. 

A suitable approximation scheme to investigate 
the quantum mechanics of this effect is provided 
by the generalized Hartree-Fock method*® worked 
out in connection with the theory of superconduc- 
tivity. The pairing forces are represented in this 
scheme by a self-consistent pairing potential, the 
strength of which is proportional to a gap param- 
eter 4. A reference system rotating with angu- 
lar frequency 2 introduces a coupling term in the 
self-consistent energy which is 2 times the angu- 
lar momentum operator L,. The competing ef- 
fect of the rotational motion and of the pairing 
finds its mathematical expression in the fact that 


21, does not commute with the pairing potential. 
Simple solvable models are being examined to 
investigate the detailed nature of the two simul- 
taneous effects. 

The main results of such a treatment can be 
seen from rather simple arguments. Since the 
Coriolis force counteracts the pairing correla- 
tions, the energy-gap parameter A will decrease 
with increasing 2. This is reflected in an in- 
crease in the effective moment of inertia g , and 
for sufficiently small 2 the effect may be de- 
scribed by a term in the rotational energy pro- 
portional to /*(J+1)?, where J is the rotational 
angular momentum. For larger values of the 
rotational frequency 2, when the rotational en- 
ergy is of the order of the binding energy due to 
the pairing correlation, the Coriolis force be- 
comes comparable to the pairing force and gives 
rise to major modifications in the correlations 
of the particles. For Q greater than a critical 
value 2.., the energy-gap parameter 4 vanishes 
and the particles perform an approximately in- 
dependent motion in the average potential. 

Though in the region 0<2<®,, one has a more 
involved intermediate coupling, the value 2, can 
independently be obtained from the simpler gap 
equation which results with a trial ground-state 
vector of the Bardeen, Cooper, Schrieffer type 
formed with eigenstates of the single-particle 
energy in the rotating system. In such a repre- 
sentation, the relevant matrix elements of the 
pairing interaction energy can be approximated 
by a constant which appears as an effective coup- 
ling constant so that, to second order* in, the 
coupling constant G of the pairing interaction is 
replaced by 


. (RIL 1k) \? 
B01 -40%{y; (=) ). ' (1) 


“p< 


The energies €, and matrix elements (k’|1, |k) 
refer to the one-particle states k of the nonro- 
tating system, and the averaging is over those 
states k which take part in the pairing interac- 
tion. The effect of 2 is to reduce the value G of 
the effective coupling constant, and 2, is the 
minimum value of 2 for which.the gap equation 
with this effective interaction leads to a vanish- 
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ing gap. This happens already for some finite 
value of G. For 2=0, the level spacing d is of 
the same order as the gap, and in the actual cal- 
culations carried out for the deformed nuclei? a 
reduction of G by about 30% is sufficient to cause 
4 to vanish. 

One can obtain an estimate of 2, by evaluating 
the averaged sum in (1), employing harmonic 
oscillator wave functions. For 4=0, 9= 9 rigid? 
and so the critical angular momentum J, is given 
by Je =Q¢ Irjgiq: One finds® in this way J, ~12 
for A~180, and /,~18 for A~238. These values 
are only slightly higher than the highest rota- 
tional states so far observed in these nuclei® 
(I=8 for A=180 and J=12 for A =238), and thus 
further study of high-angular-momentum states 
in these bands would be very interesting. One 
expects that in multiple Coulomb excitation 
studies, the rotational sequence based on the 
ground state will effectively terminate’ for / ~I,. 
States with J>J, will have an intrinsic structure 
that is practically orthogonal to the ground state. 
One would further expect that in nuclear reaction 
processes populating states with J>J,, the den- 
sity of states as a function of J will be governed 
by a Boltzmann factor in which the rotational en- 
ergy is given approximately by (h?/2 9 . id) I(I+1) 
+ constant, as is characteristic of a freé Fermi 
gas. 





*On leave from the Department of Mathematical 
Physics, University of Birmingham, Birmingham, 
England. 

‘An estimate of this energy may be obtained from the 
model employed in references 2 and 3. Thus, fora 
nucleus with an average level spacing d between the 
unperturbed (twofold degenerate) single-particle or- 
bitals, one has for the total binding energy associated 
with the pairing correlation, including both neutrons 
and protons, W=-A?/d, where 2A is the magnitude of 
the energy gap. Inserting the empirically observed 
values of A and d gives W ~-1.3 Mev for deformed nu- 
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clei, approximately independent of A. 

24. Bohr, B. R. Mottelson, and D. Pines, Phys. 
Rev. 110, 936 (1958); S. T. Belyaev, Kgl. Danske 
Videnskab. Selskab, Mat.-fys. Medd. 31, No. 11 
(1959); A. B. Migdal, Nuclear Phys. 13, 655 (1959); 
Soviet Phys. -JETP 37(10), 176 (1960); J. J. Griffin 
and M. Rich, Phys. Rev. 118, 850 (1960); S. G. Nils- 
son and O. Prior, Kgl. Danske Videnskab. Selskab, 
Mat. -fys. Medd. (to be published); V. G. Soloviev, 
Doklady Akad. Nauk (U.S.S.R.) 133, 325 (1960); and 
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3J. Bardeen, L. N. Cooper, and J. R. Schrieffer, 
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in the case of the noncommuting coupling terms which 
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87(2), 236 (1959). 

The use of a perturbation estimate at this point is 
justified since iQ, <<Ae, where Ae is a typical energy 
denominator in (1), i.e., the rotational motion implies 
only a small modification in the one-particle motion 
although it implies a major modification in the pairing 
correlation. Compare also with footnote 7. 

5One obtains essentially the same estimate for 2, 
from the condition that the difference between the ro- 
tational energies of the correlated and uncorrelated 
states should be equal to the energy gain -W due to the 
correlation. 

®¥For the rotational band of Hf!” and Hf'®, see, e.g., 
D. Strominger, J. M. Hollander, and G. T. Seaborg, 
Revs. Modern Phys. 30, 585 (1958); for U* see F. 
Stephens, R. Diamond, and I. Perlman, Phys. Rev. 
Letters 3, 435 (1959). 

"It should be emphasized that the effect considered 
here is quite distinct from the modifications in the ro- 
tational motion that occur when the rotational fre- 
quency becomes of the order of the intrinsic particle 
frequencies. These latter modifications effectively 
terminate the rotational bands when I ~1,,, where I, 
is the maximum angular momentum that can be con- 
structed out of the available single-particle configura- 
tions [see J. P. Elliott, Proc. Roy. Soc. (London) 
A245, 128 (1958)). For the heavy nuclei which exhibit 
well-defined rotational band structure, J,, ~ 40. 
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INTERPRETATION OF THE BERKELEY ANOMALY IN HIGH-ENERGY ? - d COLLISIONS* 


Arnold Tubis and Jack L. Uretsky 
Physics Department, Purdue University, Lafayette, Indiana 
(Received October 31, 1960) 


We should like to offer an alternate explanation 
of the “bump” observed in the recent Berkeley 
experiments’ on double-pion production in proton- 
deuteron collisions. It will be recalled that these 
experiments consisted of an examination of the 
momentum spectrum of the final He* or H® nu- 
cleus at a fixed laboratory angle and in a momen- 
tum range corresponding to double-pion produc- 
tion. Evidence was found for a strong peaking in 
the momentum spectrum of the final nucleus at 
a momentum corresponding to a total energy of 


about 310 Mev of the two pions in their own center- 


of-momentum system. Abashian et al. found that 
they could fit this peaking to a p-wave resonance 
in the T=1 state of the outgoing meson pair. 

We shall proceed according to the rules set up 
in reference 1; that is to say, we shall consider 
a generalization of the statistical model to allow 
for p-wave scattering between the final-state me- 
sons. We do this by writing the production matrix 
element,’ in the system where the center-of- mo- 
mentum of the two pions is at rest, in the form 


aw’) =Mq-X,, (1) 


where M is a constant, q the momentum of one of 
the pions, and X, the momentum of the final nu- 
cleus. Final-state scattering is accounted for by 
writing the total matrix element as 


g1u(q) = fea v*@,q’)w’@’). (2) 


The wave function ~*(q, q’) describes the elastic 
1-1 scattering from a state of relative momen- 
tum q’ to one of momentum q. We use the asymp- 
totic form of the wave function (in configuration 
space) and invoke a minimum “radius of inter- 
action,” R. The transition probability is then 
found to be 


PW/dpdQ, « p,7r,°E,*t[(t - 4) /t 


x (g* cotd, - R~*)* (g® cot*6, +¢%)~*. (3) 


The last equation introduces p, and E,, the labo- 
ratory momentum and energy, respectively, of 
the outgoing nucleon, ¢ is the square of the en- 
ergy of the two outgoing pions in their own rest 
system in units of the pion mass, and 6, is the 
1-a p-wave scattering phase shift for a c.m. 


energy of Vt. Let us now make the simplest possi- 


ble assumption concerning cot5,, namely that the 
scattering length approximation is valid over 
most of the energy range in question. We write, 
then, that 


q® coté, =a~’, (4) 


and choose the scattering length, a, to give the 
production cross section a maximum at about 
300-Mev total energy in the two-pion rest sys- 
tem. The scattering length turns out to be about 
2.5 meson Compton wavelengths. The corres- 
ponding momentum dependence of the outgoing 
nucleus (Tproton= 743 Mev) is shown in the solid 
curve of Fig. 1. The dashed curve, shown for 
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FIG. 1. Contribution of the T =1 two-pion production 


mode to the He® or H® momentum distribution function 

as given by Eqs. (3) and (4). The solid curve corres- 
ponds to a scattering length a of 2.5 pion Compton wave- 
lengths while the dashed curve corresponds to a zero 
scattering length. The distribution functions are evalu- 
ated for an incident proton energy of 743 Mev in the lab- 
oratory and are to be compared with Figs. 3 and 4 of 
Abashian et al.! We have used arbitrary units for the 
ordinate scale. 
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purpose of comparison, corresponds to the “p- 
wave statistical model” obtained by allowing the 
scattering length to vanish. 

It is our feeling that the present Berkeley data 
are consistent with the scattering-length model 
we have presented and that there is no necessity 
for invoking the existence of a resonance to ex- 
plain the Berkeley “bumps.” The validity of our 
approximation rests upon the assumption that the 
effective range and higher terms in the expansion 
implied by Eq. (4) are sufficiently small to war- 
rant their neglect in the energy range under con- 
sideration. We note that if this is the case, and 
if the appropriate choice of signs is made for the 
effective-range term, it is possible, within the 
framework of our model, to have a high-energy 
resonance such as the one required by Frazer 
and Fulco® although its shape will not satisfy their 
requirements. 

If the explanation we propose is found to be con- 
sistent with future experiments, it raises the 
possibility of an interesting speculation. We re- 
call that* if there is a bound p state of binding en- 
ergy Ep, and if the terms involving higher powers 
of qg? in the expression for cot6, are sufficiently 


small, then the binding energy is related to the 
scattering length by the relation (i = p =c =1) 

E B = 1/ a. 
Thus, if the scattering length in our model can 
be shown to be positive, we are led to considera- 
tion of a particle that is a two-meson bound state 
of spin and isotopic spin 1 and of a mass about 
equal to 450 electron masses.* 





*This work was supported in part by the Air Force 
Office of Scientific Research and the U. S. Atomic En- 
ergy Commission. 

‘A, Abashian, N. E. Booth, and K, M. Crowe, Phys. 
Rev. Letters 5, 258 (1960). 

*Julius S. Kovacs, Phys. Rev. 101, 397 (1956). 

5w. R. Frazer and J. R. Fulco, Phys. Rev. Letters 
2, 365 (1959). 

‘L. D. Landau and E, M. Lifshitz, Quantum Me- 
chanics (Pergamon Press, New York, 1958), pp. 401-2. 
sy. L. Uretsky and T. R. Palfrey (to be published) 

have considered the possibility of finding a two-meson 
bound s-state in photoproduction. In that paper the 
conclusion is reached that there is at present no persua- 
sive evidence contradicting the existence of a two-meson 
bound state. 





RESONANT SCATTERING OF ANTINEUTRINOS 


J. C. Barton 
Department of Physics, Northern Polytechnic, London, England 
(Received October 27, 1960) 


It has recently been suggested by Glashow’ that 
there might be a resonance in the scattering of 
antineutrinos by electrons which would lead to an 
appreciable production of 1 mesons. It was esti- 
mated that (under the hypothesis of an interme- 
diate charged boson of mass ™ 7) the antineutri- 
nos in the cosmic radiation would produce a mes- 
on flux of from 0.1 m~* day~ (for m z = nucleon 
mass) to 2 m~* day”* (for mz =K-meson mass); 
this rate should be independent of the depth be- 
low ground and should thus be distinguishable 
from mesons produced in the atmosphere. 

While studying the intensity of cosmic radiation 
at various depths underground, a recorder of 
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geometric area 0.08 m? has been operated for 21 
days at a depth of 5050 m.w.e. No cosmic-ray 
events were observed. The effective area of the 
recorder for high-energy mesons, which would 
be accompanied by knockon and pair-produced 
electrons, is at least 25% greater than the geo- 
metric area so that the meson flux is probably 
less than 0.5 m~? day~*. According to Glashow’s 
theory this result makes it unlikely that there is 
a resonance involving a virtual charged boson 
with mass equal to that of a K meson. 


‘Ss. L. Glashow, Phys. Rev. 118, 316 (1960). 
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MUONIUM FORMATION IN SEMICONDUCTORS* 


G. Fehert 
Bell Telephone Laboratories, Murray Hill, New Jersey and Columbia University, New York, New York 


and 


R. Prepost and A. M. Sachs 
Columbia University, New York, New York 
(Received November 4, 1960) 


A well-known consequence of parity noncon- 
servation in 7 - » decay’ is the longitudinal spin 
polarization of muons in meson beams emerging 
from synchrocyclotrons.? The residual polari- 
zation of positive muons before decaying in vari- 
ous solids in which they have been stopped has 
been measured by several groups.*~* In some 
substances (e.g., graphite, metals) no depolari- 
zation occurs; in others (e.g., nuclear emulsions, 
polyethylene) a partial depolarization has been 
observed; while in a few (e.g., SiO,) there is 
complete depolarization. The details of the pro- 
cesses involved are at present not understood. 
However, measurements including the quenching 
of depolarization in external longitudinal fields*"™* 
are consistent with a mechanism in which the 
muon undergoes short-lived attachments to elec- 
trons; i.e., ute~ systems called muonium are 
formed. In order to gain a better understanding 
of these processes we have studied the depolari- 
zation in semiconductors, the electronic structure 
of which is well understood, and the electronic 
concentration of which can be easily varied. In 
this Letter we report on the depolarization of 
muons in silicon and germanium, and in parti- 
cular show that by varying the concentration of 
impurities, the entire range between maximum 
and essentially zero polarization can be covered. 

The polarization of the muons is detected by 
measuring the asymmetry in the angular distri- 
bution of positrons from the muon decay. This 
is accomplished by a standard precession-type 
experiment, as illustrated in Fig. 1. Five scintil- 
lation counters were placed in the “85-Mev” ex- 
ternal positive-meson beam of the Nevis synchro- 
cyclotron. Sufficient absorber was placed between 
counters 1 and 2 to absorb the pions and allow 
the muons in the beam to stop in the target ma- 
terial placed between counters 3 and 4. Forward 
decay positrons in counter pair 23 and backward 
positrons in 45 were counted during a gate 1.2 
usec long, which was started 0.6 psec after each 
stopped muon. Alternate measurements were 
taken in one case with no applied external field 
and in another where the muon spins were pre- 
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FIG. 1. Experimental arrangement for measuring 
the depolarization of muons. Cloud-chamber magnet 
coils are used in experiments to measure the quenching 
of muon depolarization in large longitudinal fields. 
Helmholtz coils provide the transverse field in which 
the muons precess through 180 deg. 


cessed through an average angle of 180 deg ina 
40-gauss vertical field provided by a set of Helm- 
holtz coils. 

Each target was assembled from cylindrical 
samples with the desired impurity concentration. 
The effective target sizes ranged from 23in. x2+4in. 
xin. to 5in. x 5in. xlin. The targets were imbedded 
in Styrofoam holders which positioned the sam- 
ples in addition to serving as Dewars for cooling 
to liquid nitrogen temperatures. 

The experimental results are shown in Fig. 2, 
where the usual positron asymmetry parameter, 
a, is plotted against the concentration of free 
electrons and holes in silicon. The results on a 
single sample of n-type germanium and of graph- 
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FIG. 2. Experimental values of the asymmetry parameter, a, for decay positrons 


from stopped muons (a) versus free electron concentration in n-type silicon and free 
hole concentration in p-type silicon; (b) in one sample of n-type germanium (phosphorus- 
doped) at room temperature and liquid nitrogen temperature; and (c) in a graphite sample 
for which the maximum value of a=0.33 is assumed to correspond to full muon polari- 
zation. The abscissas for n-type and p-type silicon have been joined at the value of the 
intrinsic concentration for room temperature (~10'° cm’). Since the product of the 
numbers of free holes and electrons in thermal equilibrium with the lattice is constant 

at a given temperature (i.e., ~ 10” for silicon at room temperature), the entire ab- 
scissa represents an increasing free electron concentration to the right (or an increas- 


ing hole concentration to the left). 


ite are shown on the same graph. The angular 
distribution of decay positrons is of the form 
1-acosé@, where @ is the angle between the di- 
rections of the incident muons and the emitted 
positrons. In order to obtain a from the meas- 
ured peak-to-valley ratios for forward and back- 
ward positrons, corrections have been made for 
accidental counts, for positrons not originating 
in the target, and for smearing caused by the 
finite gate width used and acceptance angle sub- 
tended by the detectors. 

Before discussing the results, it may be re- 
called that embedding impurity atoms such as 
phosphorus in a material with a large dielectric 
constant €, such as silicon, creates what are 
called “shallow donors.”™* The donor impurities 
have hydrogen-like orbits of greatly expanded 
radius, an ionization energy reduced by ~«€?, and 
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markedly reduced hyperfine interactions. For 
example, phosphorus in silicon has an ionization 
energy £;=0.045 ev, and the hyperfine interaction 
of its electron is approximately 30 times smaller 
than it would be if « were unity.** Therefore, if 
muonium is formed in silicon or germanium, one 
might expect it to behave like a shallow donor, 
with reduced ionization energy and hyperfine inter- 
action. The experimental results will be discus- 
sed with this point of view in mind. 

The free electrons and holes are supplied by 
the added impurities which are almost all ionized 
at room temperature.’* The time 7g that the elec- 
tron is bound to a donor (i.e., the lifetime of the 
muonium state) can be shown to be 


T p= (ovN)* exp(E ,/kT), (1) 
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where o is the capture cross section of an ionized 
donor, v is the velocity of an electron, Nc is the 
density of states in the conduction band, E; is the 
ionization energy of the impurity, and kT is the 
thermal energy of the carriers.'® For donor con- 
centrations exceeding 10** per cm*, the donor 
wave functions overlap, the electrons are no 
longer localized, the effective 7g becomes much 
shorter than given by Eq. (1), and one approaches 
a state similar to a metal (see Fig. 16 of refer- 
ence 8). This situation is exemplified by sample 
No. 5 which like graphite shows no depolarization, 
corresponding to the lack of muonium formation. 
As the donor concentration is reduced, the elec- 
trons are able to attach themselves for a time 

Tp tothe muon. This results in the depolariza- 
tion exhibited by samples 2, 13, 8, 12, 1, and 11. 
As the donor concentration is further decreased 
and acceptors are added, the equilibrium number 
of electrons becomes too small to have a signifi- 
cant probability to form muonium.’* The partial 
depolarization observed (see samples 10, 14, 

and 6) must therefore be due to the capture of 
electrons that are not in equilibrium with the lat- 
tice but may be created, for example, by the 
muon passing through the sample. The weak de- 
pendence of a on the hole concentration is at pre- 
sent not understood. 

A germanium sample with ~10** elec/cm* 
showed only a very small degree of depolariza- 
tion (see Fig. 2). This presumably results from 
the fact that the ionization energy of donors in 
germanium is four times smaller than in silicon, 
resulting in a decrease in 7p [see Eq. (1)]. By 
going to 77°K, the relevant ratio of E;/kT equals 
that of silicon at room temperature and indeed 
an almost complete depolarization (similar to that 
for silicon) is observed (see Fig. 2). 

The polarization may be re-established by ap- 
plying a longitudinal magnetic field to the sam- 
ple.*-"4 We have performed a preliminary experi- 
ment on sample 11 and found that a field of approxi- 
mately 1000 gauss quenched half of the muon de- 
polarization. This is consistent with the idea of 
shallow donor formations if the number of (y*e~) 
attachments with the reduced hyperfine interaction 
formed in a muon lifetime is of the order of a 
hundred.’” Further depolarization experiments 
at low temperatures are in progress. 
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PHOTOPRODUCTION OF NEGATIVE PIONS FROM HYDROGEN AT FORWARD ANGLES* 


J. R. Kilner, R. E. Diebold, and R. L. Walker 
California Institute of Technology, Pasadena, California 
(Received November 3, 1960) 


The recent theoretical investigation by Drell’ 
of the photoproduction of very high energy par- 
ticles at small angles predicts a large peak in 
the differential cross section in this region. In 
particular, if one applies his results to the photo- 
production of negative pions from hydrogen one 
obtains a value for the differential cross section 
that is several times larger than the value meas- 
ured in pair production experiments at larger 
angles and lower energies. This interesting re- 
sult has prompted the experimental investigation 
reported here. 

The bremsstrahlung beam from the CalTech 
synchrotron was used to produce negative pions 
in a liquid hydrogen target predominantly through 
the process 


y+p—-a- +n +p. (1) 


This is the only process kinematically possible 
for Q less than 140 Mev, where Q is the total ki- 
netic energy in the p - 7+ center-of-momentum 
system. Above Q=140 Mev, the 7” may be pro- 
duced by 


mt+p+n° 
y+p-1 + 
mt+n+n*; (2) 


however, Chasan et al.” indicate that the cross 
section for this process is small, especially near 
the threshold. 

Since there are three or more particles in the 
final state it was necessary to perform a “syn- 
chrotron subtraction” to determine the energy, 

k, of the photon initiating the reaction. An inter- 
val of approximately 100 Mev centered at 1230 
Mev was obtained by taking the difference in nor- 
malized counting rates at peak bremsstrahlung 
energies of 1300 and 1200 Mev. 

The negative pions were analyzed with a wedge- 
type magnetic spectrometer and detected with a 
counter telescope system consisting of one Ceren- 
kov and three scintillation counters with a fifth 
“aperture-defining” counter located at the mag- 
net exit. Pole-face veto counters were used to 
eliminate particles scattering from the magnet. 
Accidental and cosmic-ray counting rates were 
negligible. Brody et al.* give a more complete 
description of the apparatus except that we used 
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the same magnet in the “high-energy” position 
as described by Vette.‘ 

Data were taken at laboratory angles ranging 
from five to thirty degrees and 1 laboratory mo- 
menta ranging from approximately 600 to 1000 
Mev/c. At laboratory angles less than about ten 
degrees appreciable numbers of electrons trav- 
ersed the system and were rejected using the fact 
that they produce showers in matter. A one-half 
inch lead converter was placed in front of each of 
the last two counters, and the “electron biases” 
on the pulses from these counters were adjusted 
so that about 95% of the electrons and only 8% of 
the pions would give a pulse above this bias. Al- 
though there were approximately three times as 
many electrons as pions at the most forward 
points measured, the largest resulting electron 
contaminations were only on the order of 20% 
because ofthe high electron detection efficiency. 

The data were further corrected for empty tar- 
get backgrounds, pion decay in flight,® 1 contami- 
nation from the 7 - » decay,°® and absorption in 
the system.® The photon beam was monitored by 
ionization chambers calibrated with an energy- 
independent Quantameter.® The error in this cali- 
bration and all other systematic errors (amount- 
ing to perhaps 5%) are not included in the graphs 
where only the statistical counting errors are 
shown. 

The experimental differential cross sections, 
o(8,w), for k=1230 Mev, are shown in Fig. 1 as 
a function of 7” center-of-momentum total energy 
(and also of Q) for several angles.” Not shown 
are four points for w=370 Mev (Q =293 Mev) with 
very large statistical errors, which gave values 
of o(@,w) consistent with other values near the 
lower energy end. Figure 2 gives the angular dis- 
tribution for the average value of the data points 
at w =464, w=505, and w=519 Mev. The peak 
value of the cross section is approximately five 
times larger than the pair production cross sec- 
tions at lower energies and backward angles.® 
However, the peak is rather narrow in both en- 
ergy and angle so that it contributes only about 
25% of the total cross section. Specifically, the 
region of the peak, 10°<@ < 55° and 420 <w < 573 
Mev in the center-of-momentum system, contri- 
butes 13 yb to the total cross section of about 50 
ub.?»® 
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FIG. 1. The differential cross section as a function 
of w (or Q) for several fixed angles. o(@,w)dQdw is 
the cross section for producing a 7~ at angle @ within 
dQ and total energy w within dw, all in the center-of- 
momentum system. The solid curves are drawn to fit 
the experimental points. The dashed curves are cal- 
culated from Drell’s model (for @=17.2 degrees) and 
for an isotropic statistical model, which are described 
in the text. Experimental resolutions have been folded 
in. 


Drell’ gives the diagram shown in Fig. 3 as the 
main contributor at forward angles for high-en- 
ergy 1 photoproduction. The differential cross 
section is given by® 


a sin?0 dQ wk - w)dw 


0(@ am Pome 
(0, w)d du 2n (1 - 8 cosé@)? 4n RS 





at++p 


(3) 


The case when more than one pion emerges from 
vertex (b) of Fig. 3 is included since this is con- 


(Q). 
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FIG. 2. The angular distribution for the average 
value of points at w=464, w=505, and w=519 Mev. 
The prediction of the Drell model averaged in the same 
manner is also shown with the experimental resolutions 
folded in. 








FIG. 3. The diagram for the photoproduction of 
negative pions which gives a large contribution at high 
energy and forward angles. Equation (3) gives the 
cross section from this diagram. 


tained in the total cross section,” o7+ 4 p(Q). 
Drell’s process may be particularly easy to recog- 
nize in hydrogen where it should show not only 

the characteristic forward angle peak, but also 

a strongly peaked energy dependence arising from 
the 3-3 resonance in o7+ + p(Q). 

The cross section calculated from Eq. (3) is 
shown in Fig. 1 for 6=17.2 degrees and in Fig. 2 
where the same average used for the data points 
is plotted. Both of these curves have the photon 
distribution and magnet resolutions folded in. 

The data show the expected effect of the pion- 
nucleon scattering resonance which has a peak 

at Q=145 Mev. Of course, any “isobar model” 
would give a similar dependence of o(@,w) on w, 
but the fact that the angular distribution decreases 
so rapidly with angle is probab!y more specific 

to Drell’s process. The distribution predicted 
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by an isotropic statistical model™ normalized to 
a total cross section*»® of 50 ub is also shown in 
Fig. 1. 

Qualitatively, our results agree rather well 
with Drell’s model, although the experimental 
cross sections are in general considerably larger 
than those calculated from Eq. (3), and do not 
show any decrease at the smallest angle meas- 
ured. One should perhaps not expect a better 
quantitative agreement since Drell’s diagram of 
Fig. 3 is not the only one leading to pion pair pro- 
duction, and other relatively small amplitudes 
may produce appreciable effects by interference 
with the Drell term. 
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RESONANCE IN THE Az SYSTEM* 


Margaret Alston, Luis W. Alvarez, Philippe Eberhard, t Myron L. Good, 
William Graziano, Harold K. Ticho, || and Stanley G. Wojcicki 
Lawrence Radiation Laboratory and Department of Physics, University of California, Berkeley, California 
(Received October 31, 1960) 


We report a study of the reaction 
K> +p+A°+nt4+07 (1) 


produced by 1.15-Bev/c K” mesons and observed 
in the Lawrence Radiation Laboratory’s 15-in. 
hydrogen bubble chamber. A preliminary report 
of these results was presented at the 1960 Roch- 
ester Conference.’ The beam was purified by 
two velocity spectrometers.? A =° hyperon ob- 
served during the run® and the preliminary cross 
sections‘ for various K~ reactions at 1.15 Bev/c 
have been reported previously. Reaction (1) was 
the first one selected for detailed study, because 
it appeared to take place with relatively large 


probability and because the event, a 2-prong inter- 


action accompanied by a V, was easily identified. 
In a volume of the chamber sufficiently restricted 
so that the scanning efficiency was near 100%, 
255 such events were found. These events were 
measured, and the track data supplied to a com- 
puter which tested each event for goodness of fit 


to various kinematic hypotheses. The possible re- 
actions, the distribution of events, and the corres- 


ponding cross sections are given in Table I. An 
event was placed in a given category of Table I if 
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the x” probability for the other hypotheses was 
<1%. It appears likely that the majority of the 
events in group (e) are also reactions of type (1). 
This belief is based on the following arguments: 
1. Since the kinematics of a Ann fit (four con- 
straints) are more overdetermined than those of 
a D°nn fit (two constraints), it is relatively easy 
for a Anz reaction to fit a 5°mm reaction, but only 


Table I. Distribution of events among different re- 
actions. 








No. of Cross 
events section 
Reaction (mb) 
(a) K7+p—-K°+p+a- 48 2.0+0.3 
(b) K~+p—(A or 2°) +n*+n7-+n" = 39 1.1+0.2 
(c) K~+p—~2°+nt+a7 27 
(d) K> +p A+at+n- 49 4.1+0.4 
(e) K7+p—(A or 2%) +nt+07 92 
Total 255 7.2 +0.5 
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very few £° configurations can fit the Ann reac- 
tions. 

2. The events of group (e) when treated as D°n7 
reactions give a x” distribution which is much 
worse than that obtained when they are treated as 
Ann reactions. 

In what follows, the 141 events of groups (d) 
and (e) are treated as examples of reaction (1). 
We estimate that 10 to 15% are actually £° events. 

The energy distribution of the two pions in the 

K~ -p barycentric system is shown in Fig. 1. If 
the cross section were dominated by phase space 
alone, the distribution of the points on the two- 
dimensional plot of Fig. 1 should be uniform. 
This is clearly not the case. On the contrary, 
both the 7+ and the 2~ distributions have peaks 
near 285 Mev, such as would be expected from 
a quasi-two-body reaction of the type 


K~+p+Y*t+4q", (2) 
the Y* having a mass spectrum peaking at ~ 1380 


Mev. If the Y* of mass 1380 Mev breaks up ac- 
cording to 


y** ~ A% at, (3) 


the pions from this breakup are expected to have 
energies ranging from 58 to 175 Mev in the K p 


rest system. Those pions from (3) are well sepa- 
rated from the pions arising from reaction (2) in 
the energy histograms. 

The isotopic spin of this excited hyperon must 
be one, since it breaks up into a A anda vz. Since 
the Y* is produced with a pion, also of isotopic . 
spin one, the reaction could proceed either in 
the J=0 or the J=1 state. Therefore the ratio of 
Y** to Y*~ will depend on the relative magnitude 
and phase of the two isotopic-spin amplitudes and 
thus could differ from unity. We observed 59 Y** 
events and 82 Y*~ events, using the criterion for 
separation that the high-momentum 7 meson is 
the pion from reaction (2). 

Figure 2 shows the distribution in mass of the 
Y* state (both Y** and Y*~) including all 141 
events, again using the higher energy pion in 
each event to calculate the Y* mass. The experi- 
mental uncertainty in the mass for each event is 
small compared to the observed width of the 
curve. The curves of Fig. 2:are discussed later. 

Figure 3 shows production angular distributions 
for Y** and Y*~ in the Kp rest system. Partial 
waves with /> 0 appear to be present, as would 
be expected since hk/m,c approximately equals 3. 
The difference between the Y** and Y*~ angular 
distributions may reflect the different superposi- 
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FIG. 2. Mass distribution for Y* and fitted curves 
for tA and mp resonances. The lower scale refers only 
to the 7A resonance. Q is the kinetic energy released 
when either isobar dissociates. The curve for the 7A 
resonances is fitted to the center eight histogram in- 
tervals of our data. The 7p curve is the fit obtained by 
Gell-Mann and Watson,’ to 7p scattering data. Both fits 
are to the formula o«X%?T?/((E -E,)?+4 r*], where I 
= 2b(a/K)3/[1 + (a/K)*). 





















































tions of the isotopic-spin zero and one amplitudes 
for the two cases. 

The following two methods were used in an ef- 
fort to determine the spin of Y*. 

(a) The angular momentum of Y* was investi- 
gated by means of an Adair analysis.® We first 
restricted ourselves to production angles with 
|cos@|20.8. For this angular range the Adair 
analysis should be valid if only S and P waves 
are present in the production process. We then 
computed 7 for each event, where 

n=P,--P,/(IP,-ll PA 1). 
Of the 29 events with |cos@|20.8, the fraction 
0.62+0.09 has |7|20.5. If the above-mentioned 
restriction on the angular interval is sufficient 
to insure the validity of the Adair analysis, this 
ratio is expected to be 0.50 for j =1/2 and 0.73 
for j=3/2. The experimental result is thus~1.3 
standard deviations from both possibilities, and 
no conclusion may be drawn from the data. Simi- 
lar results were obtained for several larger val- 
ues of the cutoff angle. Presence of D waves, 
however, cannot be excluded by the production 
angular distributions (Fig. 3). If they are pre- 
sent, indeed, then none of these choices of angle 
would be sufficiently restrictive to guarantee the 
success of the Adair analysis. 

(b) Since Y* may be polarized perpendicular to 
its plane of production, correlations can exist be- 
tween the decay angle of the Y* and the polariza- 
tion of the resulting A. Also, a net A polarization 
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can result. With our limited data, we see no 
statistically significant A polarization or angular 
correlations. 

However, one can also look for anisotropy, i.e., 
a polar-to-equatorial ratio, in the decay angle of 
Y* with respect to the normal to the plane of pro- 
duction. For spin 3/2, the distribution must be 
of the form A + Bé? by the Sachs-Eisler theorem,® 
independent of the Y* parity, where we have 


§=(Py-XPy«)*P, (IP, -xPys IP. 


and Py is the momentum of the particle a in the 
Kp barycentric system. 

Since the coefficient B is a function of the pro- 
duction angle, we want to restrict ourselves to 
that range of the solid angle where the polar-to- 
equatorial anisotropy is probably greatest along 
the normal to the production plane. For produc- 
tion angles near 0 deg and 180 deg (Adair-analy- 
sis region), one expects the polar-to-equatorial 
ratio to be most different from unity in another 
direction (namely along the direction of the beam). 
Thus the equatorial region of production angles 
is more likely to show a large anisotropy along 
the direction in question. Therefore the produc- 
tion-angle range sin@ > 0.866 was selected for 
study. We find the ratio of events with |&|>0.5 
to all events is 0.355. If the distribution is iso- 
tropic, as is required for spin 1/2, we expect 
0.500+ 0.063 for our 62 events. The result is 
thus 2.3 standard deviations from isotropy. The 
45-to-1 odds against isotropy overstate the case 
for higher spin because this is the fourth aniso- 
tropy looked for. 

Since Y* may be regarded as a hyperon isobar, 
which decays into a 7 anda A, it evidently cor- 
responds to a resonance in pion-hyperon scatter- 
ing. The mass distribution of Fig. 2 then invites 
a comparison to the cross section for pion-nu- 
cleon scattering in the 3/2 - 3/2 state. For this 
purpose a p-wave resonance formula employed 
by Gell-Mann and Watson’ for pion-nucleon scat- 
tering was fitted to our 7A data by using the eight 
central histogram intervals of Fig. 2. In fitting 
the curve, it was found that the interaction radi- 
us (a2) could be varied over a wide range without 
changing the goodness of fit appreciably, provided 
that the reduced width (b) was also changed appro- 
priately. The radius parameter was therefore 
fixed arbitrarily at h/m,c. Table Il summarizes 
our results for Y*, along with those of Gell-Mann 
and Watson for the 3-3 resonance. 

Even if Y* does turn out to be a p-wave reso- 


Table Il. Parameters for -A and 1 -p resonance 
fitted to ox xX? T?/[(E -E,)*+41%4, where Fr =2b(a/x)3/ 
[1+ (a/A)*). 





Parameter nm-p m-A 





Interaction radius a 


in units of i/myc 0.88 1 
Reduced width 6 (Mev) 58 33.4 
Resonance energy Ey (Mev) 159 129.3 
Full width at half maximum (Mev) 100 64 


Lifetime (sec) ~ 10-3 





nance, there are still many reasons why the 7-A 
resonance parameters must not be taken too liter- 
ally: (a) There is a small contamination of 2°n7 
events in our data. (b) A nonresonant background 
may be present. (c) The production matrix ele- 
ment for reaction (2) might well depend on the 
outgoing momentum, and hence distort the mass 
distribution of Y*. (d) Two thresholds for other 
possible decay modes of Y* appear within the 
mass interval covered by the resonance curve; 
i.e., the 2a mode threshold around 1330 Mev and 
the KN threshold around 1435 Mev. This must 
have some effect on the shape of the mass spec- 
trum as observed via the Az decay mode. (e) Final- 
state pion-pion interaction could disturb the spec- 
trum. (f) Even when the two resonances, Y** 

and Y*~, are well resolved in terms of intensity— 
as in our experiment—there can still be an appre- 
ciable interference between the amplitude in which 
the n+ arises from reaction (2) and the 1- from 
reaction (3) and the amplitude in which the roles 
of the two pions are reversed. 

If we bear all these uncertainties in mind, the 
resemblance to the 3-3 resonance is certainly 
remarkable (Fig. 2). The resonance energies 
when expressed in terms of barycentric kinetic 
energies differ by only 30 Mev, which is much 
less than the width of either resonance. Further- 
more, the widths are at least comparable. 

These results are strongly reminiscent of the 
concept of global symmetry which predicts two 
spin 3/2 pion-hyperon resonances, one with T=1, 
the other with T=2.* These are the hyperon | 
counterparts of the J =7=3/2 resonance of the 
pion-nucleon system. On the other hand, the pos- 
sibility that Y* is a J=1/2 resonance cannot be 
excluded on the basis of our data. The concept 
of pion-hyperon resonance in either J =1/2 or 
3/2 state has been discussed recently by several 
authors.® 
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A study of 5* 7* n° events in our experiment is 
under way at present. The results, however, are 
too incomplete for us to be able to draw any de- 
finite conclusions. 

The authors are greatly indebted to the bubble 
chamber crew under the direction of James D. 
Gow for their fine job in operating the chamber, 
especially Robert D. Watt and Glen J. Eckman 
for their invaluable help with the velocity spectro- 
meters. We also gratefully acknowledge the co- 
operation of Dr. Edward J. Lofgren and the Beva- 
tron crew, as well as the skilled work and co- 
operation of our scanning and measuring staff. 
Special thanks are due the many colleagues in our 


group who developed the PANG and KICK computer 


programs-—especially Dr. Arthur H. Rosenfeld, 
and to Dr. Frank Solmitz for many helpful dis- 
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SOME CONSIDERATIONS ON THE RECENTLY FOUND EVIDENCE FOR A 7A RESONANCE 


D. Amati and B. Vitale 
CERN, Geneva, Switzerland 


and 


A. Stanghellini 
Istituto di Fisica dell’Universita, Bologna, Italy and Istituto Nazionale di Fisica Nucleare, Sezione di Bologna, Italy 
(Received November 7, 1960) 


Some interest in the possible existence of reso- 
nant states in the pion-hyperon system has been 
raised of late by the analysis of the pion spectra 
in the reaction K"+p~A+n++z~.' The experi- 
mental data fit indeed a phenomenological picture 
where about 75% of the 7* in the previous reac- 
tion are produced together with a resonant 7A 
state, whose @ value is peaked at about 115 Mev 
and whose half-width is ~30 Mev. A detailed ana- 
lysis of such a resonant state is still missing and 
many of its characteristics (like its total and 
orbital angular momentum and the branching ra- 
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tios in decay) are as yet unknown. It is, however, 
clear that the very presence of a resonant pion- 
hyperon state is of importance, as it will help in 
a deeper understanding of the processes among 
strongly interacting particles.” 

On the other hand, the existence of this reso- 
nance can give us more confidence in the validity 
of some theoretical models for the investigation 
of the pion-hyperon interaction. This was indeed 
the case with the one-pion approximation static 
model of Chew and Low for pion-nucleon scatter- 
ing. As a matter of fact the subsequent dispersive 





i- 


ve 





VoLUME 5, NUMBER 11 


PHYSICAL REVIEW LETTERS 


DECEMBER 1, 1960 





approach made us understand in a clearer way 
the fact that just the large 3-3 resonance, by en- 
hancing the one-pion intermediate state, makes 
the other contributions relatively less important 
in low-energy pion-nucleon physics. 

An analogous model for pion-hyperon interaction 
was studied some time ago for the case of equal 
A and © parities;® it is the purpose of this Letter 
to examine again its results and predictions in 
the light of the present evidence for the 7A reso- 
nance. In our case the presence of possible KN 
intermediate states does not make the one-meson 
approximation unreliable, as the mass of the KN 
states is higher than that of the Anz states, and 
besides no enhancement or resonance is known in 
low-energy KN scattering. 

The results of our model, in which the mass 
difference between A and = is taken into account 
explicitly, depend on three parameters: the two 
normalized pion-hyperon coupling constants, fa 
and fs, and a parameter 2 which is related to the 
cutoff energy. The binding energy of A in nuclear 
matter seems to indicate that f, is of the order 
of the pion-nucleon coupling constant; nothing 
more can be said at present about the numerical 
values of f, and fs. We shall discuss in what 
follows the predictions of our model for vaiues of 


B= (fF .*-FeWAS 4? +h y?) (1) 


smaller than, let us say, ~0.3 (which means 
0.5fy? < fa? < 2f57). At the end we shall briefly 
indicate what conclusions can be obtained when 
more extreme relations between the two coupling 
constants are assumed. 

With the limitation of 6 just mentioned, our 
model predicts two resonant states for J=3/2, 
with T=1 and T=2, respectively; no resonance 
can appear either in the J=3/2, T=0 or in the 
J=1/2 states. The total resonant energies are 
given by 


B'« +2 -44-4$6A 
- 24 - $04, 


E P=m,+2+$0+20A, (2) 


where A=my - ma =80 Mev. 
It is also possible to predict the half-width of 
the two resonances, which turns out to be 


Tr, = $f" 7 0.72(1 + 0.660), 


r,= fy” q°(1+0.155), (3) 


where 
fy = HS. thes (4) 


and @g and q are the pion momenta in the 7A (T =1) 
resonance and in the 7= (T =2) resonance, respec- 
tively. ‘ 

As the T =1 resonance is present both in the 7A 
and in the 72 channel, from formulas (12) in I, 
we can also calculate the branching ratio R for 
the decay of the resonant state Y * into either a 
mA or mz state. As is expected for a narrow reso- 
nance, R turns out to be nearly independent of 
the way the resonance was prepared, which means 
that 


(Zl), _ Amz InmA) = (mAlnZ=)/{mAlmtA) (5) 


at resonance. We obtain indeed 
R=(¥* +1D)/(¥Y* +A) ~4$(0.1)/A1+5)?. (6) 


Let us now briefly discuss the results contained 
in (2), (3), and (6). We note that no 62 correc- 
tions are present in (2), which makes the correc- 
tions to the resonance energies due to 6 rather 
small. It is interesting to note that, should we 
take 2.=290 Mev (as in the pion-nucleon case, 
as global symmetry would suggest), we would get 
E, =1365 Mev, which means a Q value of 110 
Mev, in quite good agreement with the experi- 
mental value for the 7A (T=1) resonance. Or 
putting it differently, we can determine the nu- 
merical value of 2 from the experimental posi- 
tion of the T =1 resonance and notice that it turns 
out to be practically coincident with the position 
of the 3-3 resonance. 

Our model predicts a J =3/2, T =2 (x=) reso- 
nance at higher energy. The difference E,? -E,* 
does not depend on 2, and is given by 


EP -E * =24 +404 =(160 +1056) Mev. (7) 


This formula represents a clear-cut prediction 
of the theory, based only on the assumption of 
equal A-2 parities and small 6. Therefore, it 
would be extremely interesting to investigate ex- 
perimentally the existence and position of a T=2 
resonance, by studying for instance the processes 
K~+p—2" +0-+n* or n++p—2++7++K° in the 
Bev region. 

As for the widths and the value of R, we note 
first that in the limited symmetry case (i.e., A 
=6=0) one would have simply 


T,=T,=i'¢, R=3. (8) 


In the expressions (4) and (6) we can therefore 
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recognize the influence of the nonzero value of A 
(besides the g* kinematical factors) and the cor- 
rections in 6. We see that while E, and the re- 
lation (5) are practically insensitive to 6, this is 
not necessarily the case for the I. 

The value of R given by (6) gives a preference 
for the resonance to decay into the 7A state, which 
seems in agreement with the experimental results, 
even if as yet we have no reliable measurement of 
the branching ratio. 

Several results of our model, obtained with posi- 
tive A-2 relative parity and small 6, are there- 
fore either in agreement with experiment or sus- 
ceptible to a direct experimental test. A positive 
check of the predictions of the model would there- 
fore give a good argument in support of the above- 
mentioned assumptions. 

Should f,? and fs” be widely different, a new set 
of resonances could develop in our model. If for 
instance fy” ~0 (we feel safer playing with fy’, 
since the order of magnitude of f,’, as already 
noted, is roughly known), then all three isotopic 
channels in the J =3/2 states would present reso- 
nances; the T=0 and T=2 would coincide at 240 
Mev above that for T=1. The J=1/2, T=1 could 
also resonate at an energy that could be of the 
same order as that corresponding to the J = 3/2, 
T=0 or 2 resonance. It is clear that an experi- 
mental investigation of the 7=0 resonance would 
‘also be of great interest. A small nonzero value 
for fs? would, however, remove the J=1/2 reso- 
nance and, less rapidly, the J=3/2, T=0 reso- 
nance. An independent idea of the relative value 
of f,? and fs” would therefore be highly desirable. 
Perhaps a way that in the near future can be ex- 
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perimentally accessible is the polology (extra- 
polation) procedure applied to the processes 
=N+AN, UN+EN in the energy region of some 
hundreds of Mev. 

In conclusion, the agreement found between one 
of the results of our model and the experiment is 
based on the assumption of equal A -= parity and 
depends, to a lesser extent, on the assignment of 
spin 3/2 to the observed resonance (as global sym- 
metry would suggest). Should these assumptions 
turn out to be not valid, we could only repeat with 
Fontenelle*: “Je ne suis pas si convaincu de notre 
ignorance par les choses qui sont et dont la raison 
nous est inconnue, que par celles qui ne sont point 
et dont nous trouvons la raison.” 

We are grateful to Dr. J. Prentki, Dr. Y. Yama- 
guchi, and Dr. Ph. Meyer for their interest in our 
work and for useful discussions. 





'M. L. Good, Proceedings of the Tenth Annual Roch- 
ester Conference on High-Energy Nuclear Physics (to 
be published). 

24 detailed analysis of some implications of the pion- 
hyperon resonance is given by Ph. Meyer, J. Prentki, 
and Y. Yamaguchi (to be published). 

3p, Amati, A. Stanghellini, and B. Vitale, Nuovo 
cimento 13, 1143 (1959); to be referred to in the fol- 
lowing as I. A similar model was studied by M. Nauen- 
berg, Phys. Rev. Letters 2, 351 (1959). For the case 
of opposite A and = parities, analogous calculations are 
in progress at CERN. In formula (12') of I, there isa 
fortuitous omission of D,°, which is given by 


D,°=1+2(5f,*-f, Iw -2f,?-f, 1d. 


‘Fontenelle, Histoire des Oracles (Paris, 1687). 
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ABSTRACTS 


In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 








NEUTRON- PRODUCING MECHANISM IN TRANS- 
VERSE PINCHES. V. Josephson, M. H. Dazey, 
and R. F. Wuerker, Physical Research Labora- 
tory, Space Technology Laboratories, Los 
Angeles, California (Received September 28, 
1960). 


Experiments performed on a transverse pinch 
assembly show that the presence of a longitudinal 
magnetic field inside of a conducting plasma and 
in the reverse direction to that of the external 
pinching field results in the generation of the 
m=0 (sausage) instability mechanism. This in- 
stability mechanism can generate high electrical 
fields which will accelerate deuterons to energies 
sufficient for neutron production in deuterium 
plasmas. Subsequent to the blowup of the insta- 
bility, the plasma-field configuration is such that 
the accelerated deuterons can continue to circu- 
late in stable orbits until lost by neutron-produc- 
ing collisions or by diffusion out of the ends of 
the geometry. 


ELECTROMAGNETIC INTERACTION OF A 

BEAM OF CHARGED PARTICLES WITH PLAS- 
MA. Jacob Neufeld and P. H. Doyle,* Health 
Physics Division, Oak Ridge National Laboratory, 
Oak Ridge, Tennessee (Received September 30, 
1960). 


The plasma -beam instability has been studied 
by Akhiezer and Fainberg under the assumption 
that 6=0, where @ is the angle formed by the 
direction of the beam and the direction of the 
growing wave resulting from the instability. 
Under these conditions the interaction is electro- 
static, i.e., the wave is longitudinal. In this in- 
vestigation the above assumption is generalized 
so as to include the case of 6 #0 and the effect of 
electromagnetic interaction. For w~w,, where 
w, is the Langmuir frequency of the plasma, the 
interaction is electrostatic for all values of @ 


and the resulting instability which produces a 
longitudinal wave increases exponentially in ac- 
cordance with the term exp(3V3 w,”kv, cos@/8)”* 
(where w, is the Langmuir frequency of the beam). 
For a frequency range below w, the instability 

is less pronounced. However, this instability is 
significant, since the interaction is electromag- 
netic and the “growing wave” resulting from this 
interaction is characterized by an electric vector 
having both transverse and longitudinal compo- 
nents. In investigating the above instabilities an 
assumption was made that the density of the in- 
cident beam is small and the results cover all 
values of 6 except those in the immediate neigh- 
borhood of 7/2. For 6 in the neighborhood of 

1/2 the assumption is more general and the re- 
sults apply to any density of the beam. 


*Now at Michigan State University, East Lansing, 
Michigan. 


EIGENVALUES OF DENSITY MATRICES. B. C. 
Carlson and Joseph M. Keller, Institute for 
Atomic Research and Department of Physics, 
Iowa State University, Ames, Iowa (Received 
September 8, 1960). 


For a system of N identical particles in a pure 
bound state, the density matrices of positive or- 
ders p and N-p have the same nonzero eigenvalues 
with the same multiplicities. If the number of 
nonzero eigenvalues is finite, these density ma- 
trices are unitarily equivalent. 


NEUTRON DIFFRACTION STUDY OF KRYPTON 
IN THE LIQUID STATE. Glen T. Clayton* and 
LeRoy Heaton, Argonne National Laboratory, 
Argonne, Illinois (Received September 19, 1960). 


The angular distribution of 1.05A neutrons 
scattered by krypton in the liquid state at tem- 
peratures of 117°K, 133°K, 153°K, 183°K, and 
210°K is reported. By use of a sample cell fabri- 
cated from a single aluminum crystal, it was 
possible to obtain diffraction patterns up to the 
critical temperature and pressure. The scatter- 
ing curves were Fourier-transformed to obtain 
the radial atomic distribution function at each 
temperature. The number of nearest neighbors 
calculated from these distribution functions range 
from 8.5 at 117°K to 4.0 at 210°K at distances 
ranging from 4.02 A to 4.20 A, respectively. A 
comparison of radial distribution functions from 


527 








VOLUME 5, NUMBER 11 


PHYSICAL REVIEW LETTERS 





DECEMBER 1, 1960 





arbitrarily terminated data shows that the ratio 
of the last zero atomic density position to the 
first atomic density maximum position is a func- 
tion of the amount of intensity data transformed. 
It also indicates that the width of the bowl of the 
effective potential in liquid krypton is in approxi- 
mate quantitative agreement with that of the 
Lennard-Jones 12:6 potential. 


*Present address: University of Arkansas, Fayette- 
ville, Arkansas. 


NONLINEAR INTERACTION OF AN ELECTRO- 
MAGNETIC WAVE WITH A PLASMA LAYER IN 


THE PRESENCE OF A STATIC MAGNETIC FIELD. 


I. THEORY OF HARMONIC GENERATION. R. F. 
Whitmer and E. B. Barrett, Palo Alto Labora- 
tories, General Telephone and Electronics Labo- 
ratories, Palo Alto, California (Received Sep- 
tember 16, 1960). 


The theory of electromagnetic wave propagation 
through an anisotropic ionized layer, including 
the effects of the nonlinear terms in the Boltz- 
mann transport equation, is presented. The meth- 
od of solution of the nonlinear equations involves 
an expansion of all of the dependent variables in 
a Fourier series in time. The differential equa- 
tions describing wave propagation are then solved, 
for each frequency in the series, for plane wave 
propagation, including all of the reflections within 
the plasma layer. A solution in closed form has 
been obtained, under small signal conditions, for 
the field at the hth harmonic in the Fourier series. 
A discussion of the properties of the wave at the 
second harmonic frequency as a function of the 
dc magnetic field strength, the electron density, 
the electron-neutral particle collision frequency, 
the field strength of the incident wave, and the 
thickness of the plasma layer is given. 


THEORY OF ELECTRON-PHONON INTERAC- 
TIONS. G. D. Whitfield,* Department of Physics, 
Columbia University, New York, New York, and 
International Business Machines Watson Labora- 
tory, New York, New York (Received August 30, 
1960). 


The theory of the interaction of electrons and 
acoustic phonons in nonpolar crystals has been 
formulated in terms of a new set of basis states, 
whose wave functions are essentially Bloch func- 
tions that deform with the lattice. The major 
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part of the interaction may then be calculated in 
terms of the strain tensor rather than the dis- 
placement of the lattice. A result of the theory 
is a generalization of the deformation potential 
theorem. 


*Now at Department of Physics, University of Illinois 
Urbana, Illinois. 


PHYSICAL PROPERTIES OF SEVERAL II- V 
SEMICONDUCTORS. W. J. Turner, A. S. Fisch- 
ler, and W. E. Reese, Research Laboratory, 
International Business Machines Corporation, 
Poughkeepsie, New York (Received September 
26, 1960). 


The physical properties of single crystals of 
the noncubic II - V semiconducting compounds 
Zn,As,, ZnAs,, ZnSb, Cd,As,, CdAs,, and CdSb 
have been investigated. Energy gaps in these 
materials vary from approximately 0.13 to 1.0 
ev. Mobilities at 297°K range from 10 cm?/volt- 
sec to 15000 cm?/volt-sec and increase at low 
temperature. Resistivity and mobility anisotropy 
have been investigated in detail for CdAs,, Ex- 
cept for the A, B.V compounds, high optical 
transmission has been observed from the intrinsic 
edge to approximately 30 microns. 


ROLE OF DOUBLE EXCHANGE IN THE MAG- 
NETIC STRUCTURE OF Li,Mnj.,Se. R. R. 
Heikes, Westinghouse Research Laboratories, 
Pittsburgh, Pennsylvania, T. R. McGuire and 
R. J. Happel, U. S. Naval Ordnance Laboratory, 
White Oak, Silver Spring, Maryland (Received 
August 22, 1960). 


The details of the magnetic behavior of the 
Li,Mn, _,Se system are attributed to the double- 
exchange interaction. At low temperatures, the 
hole which is introduced by the Li* is loosely 
bound to the Lit* itself. In the region of the Li 
ion double exchange causes local distortions of 
the spin system which we refer to as clusters. 
As the Li concentration is increased (x = 0.07) 
the clusters overlap sufficiently so that a mag- 
netic field will induce an appreciable magnetic 
moment (0.54). At temperatures below 45°K a 
canted spin ordering is suggested as the magnetic 
model for the x = 0.07 composition. Finally, at 
x =0.10 we find that spontaneous magnetization 
develops below 110°K. As the temperature is 
lowered through 70°K the spontaneous moment 
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disappears and antiferromagnetism is found. It 

is not inconsistent with the data for x = 0.10 that 
this antiferromagnetic state is a canted spin sys- 
tem with very small canting angle and therefore 
small magnetic moment. The theory of de Gennes 
is used in our discussion of the magnetic model. 


MAGNETIC STRUCTURE TRANSITIONS IN 
Li,Mn;.,Se. S. J. Pickart, U. S. Naval Ordnance 
Laboratory, White Oak, Silver Spring, Maryland, 
and Brookhaven National Laboratory, Upton, New 
York, R. Nathans, Pennsylvania State University, 
University Park, Pennsylvania, and Brookhaven 
National Laboratory, Upton, New York, and G. 
Shirane, Westinghouse Research Laboratories, 
Pittsburgh, Pennsylvania (Received August 22, 
1960). 


The magnetic structures occurring in lithium- 
substituted manganese selenide (Li, Mn, -,Se) 
have been examined by low-temperature powder 
neutron diffraction measurements. The composi- 
tion with x = 0.05 retains the fcc ordering of the 
second kind found in MnSe, the transition tem- 
perature being lowered to 83°K. For x =0.07 the 
same type of ordering sets in at 73°K, but the 
spin direction changes abruptly as the tempera- 
ture is lowered through 45°K; furthermore, the 
superlattice intensities decrease when an external 
magnetic field is applied along the scattering 
vector. At x=0.10, the spontaneous moment ob- 
served at 77°K by magnetization measurements 
is shown to be ferromagnetic, again by means of 
an external field, and a transition is found at 71°K 
from ferromagnetism to antiferromagnetism with 
the third kind of ordering. The results are dis- 
cussed with relation to models containing canted 
spins and multiple antiferromagnetic axes. 


NUCLEAR MAGNETIC RESONANCE IN SUPER- 
CONDUCTING TIN. G. M. Androes and W. D. 
Knight, Department of Physics, University of 
California, Berkeley, California (Received 
September 16, 1960). 


The nuclear magnetic resonance is investigated 
in small particles (~100A diam) of 8 tin between 
1.5 and 4.2°K, and in magnetic fields between 
1.2 and 8.8 kilogauss. The critical temperature 
and critical field are 3.71°K and 25 kilogauss, 
respectively. The effective penetration depth 
for the superconducting particles is estimated 


to be 1500 A. The resonance linewidth is 0.34% 
of the magnetic field, and it is independent of 
temperature. With respect to a tin, the nuclear 
magnetic resonance shift for 8 tin is 0.77% in 
the normal state; it approaches 0.59 % in the 
superconductor as T+0. (The largest known 
chemical shift is only 0.17%.) The variation 
with magnetic field is less than 0.03%. One may 
conclude that the electronic spin susceptibility 
in the superconducting particles at absolute zero 
is approximately three quarters of the normal 
value. The result for 1000A particles, though 
less accurate, is substantially the same. 


ENERGY DISSIPATION AND SECONDARY ELEC- 
TRON EMISSION IN SOLIDS. H. Kanter, Westing- 
house Research Laboratories, Pittsburgh, Penn- 

sylvania (Received September 23, 1960). 


Experimental evidence is presented for the 
proportionality between secondary electron yield 
and the energy dissipated by electrons near the 
surface of a solid. Using measurements of the 
energy carried away by electrons transmitted 
and reflected from thin foils of aluminum and 
carbon, the energy dissipated in an incremental 
layer at the exit surface was obtained. Simultane- 
ous measurements of the secondary electron yield 
showed a close proportionality between the num- 
ber of secondaries produced and the energy dis- 
sipation density near the surface independent of 
the incident electron energy between 1 and 10 kev. 
By subtracting the contribution of the back-scat- 
tered electrons to the yield at the front surface 
of a thick aluminum target, the yield of second- 
aries was found to be proportional to the rate of 
energy loss calculated from the Bohr-Bethe theory 
over the energy range investigated. 


CONTRIBUTION OF BACK-SCATTERED ELEC- 
TRONS TO SECONDARY ELECTRON FORMA- 
TION. H. Kanter, Westinghouse Research Labo- 
ratories, Pittsburgh, Pennsylvania (Received 
September 23, 1960). 


It is shown experimentally that back-scattered 
electrons emitted from solids under electron 
bombardment contribute significantly to the ob- 
served secondary yield, even for the case of low 
back-scattering coefficients. Thus, it was found 
that in Al with a back-scattering coefficient of 
only 0.14, about 40% of all secondaries are pro- 
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duced by back-scattered electrons for initial en- 
ergies from several kev to several tens of kev. 
The large contribution of back-scattered elec- 
trons to secondary formation even for low-atomic- 
number materials agrees approximately with what 
one would expect from the larger rate of energy 
loss and the greater path lengths of the back-scat- 
tered electrons in the secondary-electron escape 


region compared to that of the incoming primaries. 


EFFECT OF PRESSURE ON THE MOBILITY OF 
INTERSTITIAL OXYGEN AND NITROGEN IN 
VANADIUM. G. W. Tichelaar,* R. V. Coleman, 
and D. Lazarus, Department of Physics, Uni- 
versity of Illinois, Urbana, Illinois (Received 
September 13, 1960). 


Measurements of stress relaxation as a func- 
tion of hydrostatic pressure up to 9000 kg/cm? 
have been made on a vanadium sample containing 
approximately 0.1 atomic percent dissolved oxy- 
gen and 0.2 atomic percent dissolved nitrogen. 

In the temperature range 83.0° to 98.0°C, the re- 
laxation time due to dissolved oxygen is found to 
increase exponentially with pressure, the value 
at 9000 kg/cm? being about 1.7 times the value 

at 1 kg/cm?. The pressure dependence of the 
stress relaxation can be interpreted in terms of 
an activation volume of 1.7 cm*/mole which is 
about equal to the molar volume of the diffusing 
oxygen atoms. The relaxation times due to the 
dissolved nitrogen have been measured at 156.8°C 
and 163.0°C. The values at 9000 kg/cm? are about 
1.3 times the values at 1 kg/cm’, the activation 
volume being about 1.1 cm*/mole. 

*Present address: Philips Research Laboratories, 
Eindhoven, The Netherlands. 


IMPURITY CONDUCTION IN SILICON. R. K. Ray 
and H. Y. Fan, Purdue University, Lafayette, 
Indiana (Received September 23, 1960). 


Impurity conduction at low temperature was 
investigated for various p- and n-type silicon 
samples. Emphasis was placed on the study of 
samples of low impurity concentration where the 
conduction is attributed to charge exchange be- 
tween impurity centers which are partially ion- 
ized by some compensating impurity. A new 
method was used to determine the compensation. 
Donors were added to p-type samples by heat 
treatment. From changes in the room tempera- 
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ture resistivity and the Hall coefficient in the 
deionization range, the added compensation and 
the original compensation were determined. The 
measurement of various samples gave the de- 
pendence of the activation energy of conduction 
on impurity concentration and degree of com- 
pensation. The activation energy was much 
larger for the Ga- and Al-doped samples than 
for the B-doped samples of comparable impurity 
concentration. However, it was found that the 
high activation energies may be the result of ion- 
pairing between gallium or aluminum atoms and 
the compensating impurity in the sample. The 
conductivities of the various samples may be 
correlated by an expression of conductivity which 
involves the impurity concentration and the radius 
of impurity wave function. The results are dis- 
cussed in the light of current theories of im- 
purity conduction. 


THEORY OF THE MAGNETIC ANISOTROPY IN 
KMnF,. J. J. Pearson, Department of Physics, 
University of Pittsburgh, Pittsburgh, Pennsyl- 
vania (Received September 23, 1960). 


A theoretical calculation is made of the mag- 
netic anisotropy in the cubic perovskite struc- 
ture of KMnF, at room temperature and in its 
distorted structures at lower temperatures. 
These distortions are of two types: first, a 
small tetragonal distortion of the entire crystal; 
and then, below the antiferromagnetic Néel 
point, a distortion of the octahedron of fluorine 
atoms surrounding each manganese. The cubic 
anisotropy is obtained from a general spin-wave 
calculation of the zero-point dipole-dipole en- 
ergy in a cubic antiferromagnet. The result is 
found to be the same as that for the ferromag- 
netic case. The anisotropy from the tetragonal 
distortion is obtained from the change in the 
classical Lorentz factors. In calculating the ef- 
fect of the fluorine distortion a generalization is 
introduced of Kondo’s method for obtaining the 
anisotropic effective spin Hamiltonian produced 
by overlap and electron transfer between an Mn*t 
ion and its nonmagnetic neighbors. In its pre- 
sent form the method permits the ready calcula- 
tion of this anisotropy for any symmetry and 
number of neighbors. Comparison with the 
microwave resonance and torque measurements 
of Portis, Teaney, and Heeger reveals the last 
effect to be the most important and confirms the 
form of the spin Hamiltonian found here and its 
approximate magnitude. 
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EXCHANGE INTEGRAL IN COBALT FROM SPIN 
WAVE RESONANCE. P. E. Tannenwald and 

R. Weber, Lincoln Laboratory, Massachusetts 
Institute of Technology, Lexington, Massachu- 
setts (Received September 28, 1960; revised 
manuscript received October 17, 1960). 


The exchange constant A and exchange integral 
J, and their temperature dependences, have been 
measured in cobalt metal films by spin wave reso- 
nance. At room temperature A = 1.30x10-* erg/ 
cm® and J=155k. J is temperature independent 
between 4°K and 295°K. Comparison is made with 
recent data obtained by other experimental meth- 
ods. 


DECAY OF EXCESS CARRIERS IN SEMICON- 
DUCTORS. I. K. C. Nomura and J. S. Blakemore, 
Honeywell Research Center, Hopkins, Minnesota 
(Received August 22, 1960). 


A physical interpretation is given of the non- 
linear differential equations which govern the 
decay of excess carrier populations through re- 
combination centers. No restrictions are placed 
on the magnitudes of the excess carrier densities 
or the center density. Criteria for trapping are 
presented; with semiconductors for which the 
trapping level lies in the opposite half of the in- 
trinsic gap from the Fermi level it is shown that 
trapping can be described as being of either a 
temporary or permanent nature. The variety of 
possible modes of decay are illustrated with the 
aid of numerical solutions and approximate anal- 
ytic solutions. 


EXCESS TUNNEL CURRENT IN SILICON ESAKI 
JUNCTIONS. A. G. Chynoweth, W. L. Feldmann, 
and R. A. Logan, Bell Telephone Laboratories, 
Murray Hill, New Jersey (Received June 6, 1960; 
revised manuscript received October 20, 1960). 


At low forward biases, a high current flows in 
Esaki junctions due to band-to-band tunnelling. 
At sufficiently high biases the current flows by 
normal forward injection. Between these two 
bias ranges, the current is unexpectedly high 
and has been called the excess current. A com- 
prehensive experimental study has been made of 
this excess current in silicon junctions. It is 
shown that the properties of the excess current 
observed so far can be accounted for by a mech- 
anism originally suggested by Yajima and Esaki, 


in which carriers tunnel by way of energy states 
within the forbidden gap. Based on this model, 
the following expression for the excess current, 
I,, is proposed: 


i g Va a 
I D. exp{ (a, Wie /2)|e eV, +0.6e(V,+V,)I, 


where D, is the density of states in the forbidden 
gap at an energy related to the forward bias, V,, 
and the Fermi energies on the m and p sides are 
Vy and Vp, respectively, e is the electron charge, 
€ is the energy gap, W, is the junction width con- 
stant, and a, is a constant containing a reduced 
effective mass, m,. This formula describes the 
observed dependence of I, on (i) Dy, observed by 
introducing states associated with electron bom- 
bardment, (ii) €, studied by the temperature 
variation of the diode characteristics, (iii) Vx, 
verified from semilogarithmic plots of the for- 
ward characteristics, and (iv) W,, tested by using 
junctions of different widths. From these experi- 
ments, m, =0.3m, to within a factor of 2. 

The origins of the states in the band gap are 
not known for certain though they are most likely 
the band edge tails inherent to heavily doped 
semiconductors. It is probable that the tunnelling 
via local states model for the excess current in 
silicon is applicable to excess currents in other 
materials. 


THERMODYNAMIC GREEN’S FUNCTION METH- 
ODS IN NEUTRON SCATTERING BY CRYSTALS. 
Gordon Baym,* Lyman Laboratory of Physics, 
Harvard University, Cambridge Massachusetts 
(Received August 5, 1960; revised manuscript 
received October 27, 1960). 


Formulas are derived for the transition proba- 
bilities per unit time for both inelastic coherent 
scattering of neutrons by crystals, and resonant 
emission of photons and neutrons by nuclei bound 
in crystals, without making the assumption that 
the crystal is harmonic. In deriving these transi- 
tion probabilities, the analytic structure of thermo- 
dynamic correlation or Green’s functions, con- 
sidered as functions of complex temperatures 
and times, is developed and used. In particular 
a spectral form is found for the phonon Green’s 
function, Only one assumption is made about the 
crystal, namely that the displacement of the nu- 
clei due to the forces exerted by the neutron in 
scattering are linear functions of these forces. 
This leads to an evaluation of the transition pro- 
babilities in terms of the exact thermodynamic 
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displacement autocorrelation function. This eval- 
uation obeys the detailed balancing condition, and 
Placzek’s sum rule. A consequence of this eval- 
uation is that the widths of the “one-phonon” peaks 
in the neutron scattering are exactly equal to the 
widths of the corresponding phonon states of the 
crystal. 

*Present address: Institute for Theoretical Physics, 
Copenhagen, Denmark. 


KNIGHT SHIFT IN POTASSIUM. F. J. Milford 
and W. B. Gager, Battelle Memorial Institute, 
Columbus, Ohio (Received December 7, 1959; 
revised manuscript received September 19, 1960). 


The K** nuclear magnetic resonance has been 
observed in metallic potassium and in aqueous 
solutions of KNO, and K,Co(CN),. From these 
observations a Knight shift of (0.248+ 0.005) % 
has been obtained. Using Pines’ value 0.60x10"° 
for the spin susceptibility gives Py= 0.95 a,~°. 
This is compared with the theoretical values 
0.786 a,~* obtained by Callaway and 0.909 a,~* 
obtained by a quantum defect calculation in this 


paper. 


ELECTRON EFFECTIVE MASSES OF InAs AND 
GaAs AS A FUNCTION OF TEMPERATURE AND 
DOPING. Manuel Cardona, Laboratories RCA 
Ltd., Zurich, Switzerland (Received September 
19, 1960). 


The electron effective masses of several GaAs 
and InAs samples at room and liquid nitrogen 
temperatures have been determined from Fara- 
day rotation and infrared reflectivity measure- 
ments. An increase in effective mass with in- 
creasing carrier concentration has been found in 
both materials. This increase can be quantita- 
tively interpreted in InAs in terms of the non- 
parabolic nature of the conduction band. In GaAs 
the increase in effective mass with doping sug- 
gests the existence of another set of conduction 
band minima above the lowest (000) minimum. 
The measured temperature variation of the ef- 


fective mass can be attributed to two mechanisms: 


the increase in effective mass produced by the 
spread in the Fermi distribution because of the 
nonparabolic shape of the band, and the variation 
in the band structure produced by the thermal 
expansion of the lattice. 

The Faraday rotation due to the interband tran- 
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sitions has been measured in GaAs and InAs. 
This rotation is counterclockwise along the direc- 
tion of the radiation and the magnetic field for 
GaAs and clockwise for InAs. This effect is 
compared with the corresponding effect in other 
semiconductors. 


DRIFT VELOCITIES OF SLOW ELECTRONS IN 
HELIUM, NEON, ARGON, HYDROGEN, AND 
NITROGEN. J. L. Pack and A. V. Phelps, 
Westinghouse Research Laboratories, Pittsburgh, 
Pennsylvania (Received August 29, 1960). 


The drift velocities of electrons in helium, 
neon, argon, hydrogen, and nitrogen have been 
measured for E /p values between 10™ and 10 
volts/cm-mm Hg at temperatures between 77°K 
and 373°K. The data were obtained from meas- 
urements of electron transit time in an improved 
version of the double-shutter tube developed by 
Bradbury and Nielsen. By applying sufficiently 
small voltage pulses to the control grids, it was 
possible to eliminate end effects present in pre- 
vious experiments. Values of the momentum 
transfer cross sections for electrons with ener- 
gies between about 0.003 and 0.05 ev are obtained 
which are consistent with the measured drift 
velocities for thermal electrons in helium, argon, 
hydrogen, and nitrogen. The derived momentum 
transfer cross section for electrons in helium 
is found to be independent of electron energy and 
equal to 5.3x107* cm?. The momentum transfer 
cross sections for argon, hydrogen, and nitro- 
gen vary with electron energy. 


APPLICATION OF THE METHOD OF POLARIZED 
ORBITALS TO THE SCATTERING OF ELEC- 
TRONS FROM HYDROGEN. A. Temkin* and J. C. 
Lamkin, National Bureau of Standards, Washing- 
ton, D. C. (Received September 14, 1960). 


The s-, p-, and d-wave scattering of slow 
electrons from atomic hydrogen is calculated by 
the method of polarized orbitals. Utilization of 
a transformation of Omidvar avoids the iterative 
procedure of solving the associated integro- 
differential equations. The s-wave scattering is 
smaller than that given by the exchange approxi- 
mation, and the scattering lengths are within the 
upper bounds found by Rosenberg, Spruch, and 
O’Malley. The d-wave phase shifts are too small 
to explain a resonance in the total cross section. 
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However, they are much larger than those of the 
exchange or Born approximation, and they give 
considerable structure to the differential cross- 
section curves. The p-wave phase shifts are not 
much increased by the polarization effects. 


*Present address: National Aeronautics and Space 
Administration, Washington, D. C. 


COLLISIONS OF ELECTRONS WITH HYDROGEN 
ATOMS. VI. ANGULAR DISTRIBUTION IN 
ELASTIC SCATTERING. H. B. Gilbody, R. F. 
Stebbings, and Wade L. Fite, John Jay Hopkins 
Laboratory for Pure and Applied Science, Gen- 
eral Atomic, San Diego, California (Received 
September 28, 1960). 


The angular distribution of electrons scattered 
elastically by hydrogen atoms has been deter- 
mined for electron energies below 10 ev. The 
elastically scattered electrons arising from the 
interaction of crossed electron and modulated 
hydrogen-atom beams were examined over an 
angular range extending from 30° to 120°. The 
results are discussed with reference to other 


recent experimental and theoretical developments. 


NUCLEAR MAGNETIC INTERACTIONS IN HY- 
DROGEN FLUORIDE. Milton R. Baker, H. Mark 
Nelson,* John A. Leavitt,t and Norman F. Ram- 
sey, Harvard University, Cambridge, Massa- 
chusetts (Received September 16, 1960). 


The radio-frequency spectra corresponding to 
the reorientation of the proton and fluorine nu- 
clear magnetic moments in hydrogen fluoride 
(HF) have been observed in fields of 900, 1800, 
and 3600 gauss by means of the molecular beam 
resonance method. Details are presented on the 
design and construction of the new molecular 
beam apparatus and electron-bombardment de- 
tector used in the experiment. The theory of the 
HF strong-field energy levels is outlined, and 
the expected proton and fluorine transitions de- 
rived for J=0, 1, and 2 are tabulated. From 
the observed resonance shapes, one can deduce 
the magnitudes of the spin- rotational interac- 
tions of the proton and fluorine nuclei, and their 
spin-spin interaction. These are: lcp| =71+3 
ke/sec, |cp| =305+2 kc/sec, and d, =57+2 kc/ 
sec. The correctness of these parameters was 
checked by the good agreement between the ex- 
perimental curves and resonance shapes pre- 


dicted by Univac programs using these values. 
The observed fluorine spin-rotational interac- 
tion constant is the largest yet observed and 
corresponds to a rotational magnetic field at the 
nucleus of 76 gauss per unit rotational quantum 
number. The implications of the large spin- — 
rotational interaction for relaxation processes in 
nuclear magnetic resonance experiments are 
discussed. 


*Present address: Brigham Young University, Provo, 
Utah. 

tPresent address: University of Arizona, Tucson, 
Arizona. 


PERTURBATION TREATMENT OF HARTREE- 
FOCK EQUATIONS. Jan Linderberg,* Department 
of Physics, University of Florida, Gainesville, 
Florida (Received August 29, 1960). 


For atomic configurations (1s?), (1s?2s), and 
(1s?2s?) perturbation theory has been used to ob- 
tain (1/Z) expansions for the Hartree-Fock en- 
ergies. 

*On leave from Quantum Chemistry Group, Uppsala 
University, Uppsala, Sweden. 


INVERSE PHOTONUCLEAR REACTIONS 

N"*(p, y)O"® AND N?°(p, y)O*® IN THE REGION 

OF THE GIANT RESONANCE. S. G. Cohen,* 

P. S. Fisher,f and E. K. Warburton, Palmer 
Physical Laboratory, Princeton University, 
Princeton, New Jersey (Received September 28, 
1960). 


The 90° yield of y rays to the O** ground state 
from the N"*(p, y)O"* reaction has been measured 
for proton energies between 12 and 19.5 Mev 
covering the region of excitation in O** between 
18 and 25 Mev. The excitation curve is quite flat 
(with do/dQ at 90°=16 ub/4z7 sr), and shows little 
evidence of the giant resonance. The results for 
O'* are compared to those for N’® obtained by 
Jacobs and Stephens by means of the N**(y,p)C** 
reaction. The 90° yield of y rays to the O** ground 
state from the N**(p, y)O'* reaction has been 
measured for proton energies between 10 and 15 
Mev, corresponding to O** excitation energies 
between 21 and 26 Mev. The excitation curve 
shows two large resonances peaked at 21.8 and 
24.7 Mev, with integral total cross sections of 
about 0.27 Mev-mb each if no background is 
assumed. The O** results are compared to theo- 
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retical calculations of Elliott and Flowers and of 
others. 


*On leave from The Hebrew University, Jerusalem, 


Israel. 
tNow at Clarendon Laboratory, Oxford, England. 


HALF-LIFE OF O™. D. L. Hendrie and J. B. 
Gerhart, Department of Physics, University of 
Washington, Seattle, Washington (Received 
September 26, 1960). 


The half-life of O'* has been remeasured and 
found to be (70.914 0.04) seconds. The corres- 
ponding partial half-life for the 0* ~ 0* beta tran- 
sition is found to be (71.34+ 0.08) seconds. These 


results are about 1.6% less than reported earlier. 


The new half-life is combined with the O* end- 
point energy as recently determined by Bardin 
et al. to re-evaluate the Fermi coupling constant 
for beta decay. From this, the mean life of the 
uu meson is calculated. The predicted mean life 
is greater than the observed value by (3.3+ 0.4)% 
if the radiative corrections of Kinoshita and 
Sirlin are used, and by (1.4+ 0.4)% if the radia- 
tive corrections of Durand et al. are used. 


RELATIVE YIELDS OF NEUTRON GROUPS 
FROM THE Li"(p,n)Be’, Be’* REACTIONS. 

P. R. Bevington, W. W. Rolland,* and H. W. 
Lewis,? Duke University, Durham, North Caro- 
lina (Received September 23, 1960). 


The relative yields of the two groups of neu- 
trons from the Li’(p,n)Be’, Be’* reaction, lead- 
ing to the ground state and the 430-kev state of 
Be’, have been measured with a time-of-flight 
system, using pre-acceleration pulsing of the 
accelerator beam. Data were taken at 30° inter- 
vals between 0° and 150° for proton energies 
high enough to produce (p,n’) neutrons above 
the detection threshold (300 kev). Since the 
yield of the (p,m) reaction is more highly peaked 
in the forward direction than that of the (p, 7’) 
reaction, the ratio of the (p,n’) to (p,m) inten- 
sities grows with increasing angle, severely 
limiting the usefulness of the (p,m) reaction as 
a neutron source above the (p,n’) threshold at 
back angles. Absolute differential and total 
cross sections for both groups have been calcu- 
lated from the data. A comparison with theory 
for total cross sections and angular distributions 
suggests the existence of three previously un- 
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identified levels in Be®. One level, with J” =1-, 
near the threshold for the (~,m’) reaction, is 
responsible for the fast rise in the (p,n’) total 
cross section near threshold. A second level, 
corresponding to an incident proton energy of 
about 3.0 Mev, does not contribute significantly 
to the yield of the (p,m’) reaction; the data are 
consistent with an assignment of J” =1* anda 
total width of 1 Mev. The bulk of the total cross- 
section curve for the (p,n’) reaction has been 
fitted by assuming a 1* level corresponding to 
an incident proton energy of 3.5 Mev, with y,” 
*Yp"s Yn? 5%", and Ypr <¥p" 


*Now at Kings College, Bristol, Tennessee. 
tCurrently on leave at American University of Bei- 
rut, Beirut, Lebanon. 


ANGULAR DISTRIBUTIONS FOR C*"(d,n)N"*. 
Richard Zdanis,* George E. Owen, and L. Madan- 
sky, Department of Physics, The Johns Hopkins 
University, Baltimore, Maryland (Received 
September 2, 1960). 


Angular distributions of neutrons corresponding 
to transitions to the first (2.31-Mev) and second 
(3.95-Mev) excited states in the reaction 
C'5(d,n)N™ have been measured at an incident 
deuteron energy of 1.3 Mev. A time-of-flight 
technique utilizing the associated gamma ray as 
a trigger was employed to separate the various 
neutron groups. Because the data, particularly 
for the first excited state transition, exhibited 
a relatively large intensity in the backward di- 
rection, a simple deuteron stripping process is 
not sufficient to account for the results. An ana- 
lysis which included the heavy-particle stripping 
mode was consistent with the data. 


*Now at the Department of Physics, Princeton Uni- 
versity, Princeton, New Jersey. 


RADIATIONS FROM HIGH- ENERGY POSITRONS 
INCIDENT ON A BERYLLIUM TARGET. C. P. 
Jupiter, N. E. Hansen, R. E. Shafer, and S. C. 
Fultz, Lawrence Radiation Laboratory, Univer- 
sity of California, Livermore, California (Re- 
ceived September 23, 1960). 


The energy spread and yields of nearly mono- 
energetic photons of various energies produced 
by the annihilation in flight of relativistic posi- 
trons have been experimentally determined 
using a 6-in. long by 5-in. diam Nal(Ti) crystal 
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spectrometer. Photon lines with an energy spread 
of a few percent are reported. Yields of positrons 
and monoenergetic photons in the energy range 
from 2 Mev to 14 Mev have been measured. The 
bremsstrahlung spectra from 8.5- Mev positrons 
and electrons were compared and the yields and 
spectral distributions were found to agree within 
the experimental error. 


LITHIUM-INDUCED REACTION YIELDS BELOW 
4 Mev. Edwin Norbeck,* School of Physics, Uni- 
versity of Minnesota, Minneapolis, Minnesota 
(Received September 19, 1960). 


Thick-target yields of the following reactions 
were measured by counting the beta-active pro- 
ducts: Li’(Li®,2n)C™“, Be®(Li®, 2n)N'*, C’#(Li®, 
n)F"", c?(Li’, 2) 7. N**(Li®, He*)o**, N"*(Li®, 
d)F'®, N**(Li’, )F’®, O'*(Li®,2)Na”™, O7*(Li®, 
He*)F’®, F’°(Li®, Li’) F*, Na?*(Li®, Li®)Na™. The 
reactions F’*(Li®, 2p)Ne”* and Na?°(Li’, Li®)Na™* 
had too small a yield to permit accurate meas- 
urement. All of the yield curves show a very 
rapid but smooth increase of yield with energy. 
Some general rules are given for estimating the 
yield to be expected for any positive Q lithium 
beam reaction in the energy range under con- 
sideration. 


*Now at State University of Iowa, Iowa City, Iowa. 


PROTON GROUPS FROM THE F°%(a, p)Ne*” RE- 
ACTION AND THE Ca*(a,p)Sc*® REACTION. 

H. J. Martin,* M. B. Sampson, and D. W. Miller, 
Indiana University, Bloomington, Indiana (Re- 
ceived September 26, 1960). 


A CaF, target was bombarded with 21.9-Mev 
alpha particles and the energies of the outgoing 
protons were measured with a magnetic spectro- 
meter. Proton groups leading to energy levels 
in Ne*”* at excitation energies of 0, 1.28, 3.37, 
4.52, 5.18, 5.67, 6.41, 6.88, and 7.49 Mev 
were seen. The ground-state Q value for the 
Ca*°(a, p)Sc** reaction was found to be -3.47 
+0.030 Mev. Some information about absolute 
cross sections and angular distributions for the 
fluorine and calcium reactions is also presented. 


*Now on leave at Brookhaven National Laboratory, 
Upton, New York. 


PHOTONUCLEAR STUDIES WITH MONOENER- 
GETIC GAMMA RAYS FROM THERMAL NEU- 
TRON CAPTURE. Robert E. Welsh* and D. J. 
Donahue, Department of Physics, Pennsylvania 
State University, University Park, Pennsylvania 
(Received September 28, 1960). . 


Monoenergetic y rays produced when neutrons 
are captured in various materials have been used 
to study the following reactions: Ta’®"(y,n)Ta!®”” 
(8.15 hr), Au’®"(y,”)Au’® (5.6 days), Ho*®*(,,)Ho'™ 
(34 min), Ag’*"(y,m)Ag'®* (24 min), and Nb®*(y,n)Nb™ 
(10 days). Cross sections for these reactions 
have been obtained at several discrete energies 
between 7.5 and 10.8 Mev. Estimates were also 
made of the thresholds of these reactions. Two 
of the thresholds so obtained Ta’®*(y,n)Ta!®™, 
E;=17.60+ 0.08 Mev, and Nb*(y,)Nb™, E;=8.99 
+ 0.04 Mev, have precision comparable with those 
of previous measurements of the same quantities. 


*Present address: Carnegie Institute of Technology, 
Pittsburgh, Pennsylvania. 


THREE-BODY NUCLEAR PROBLEM WITH RE- 
PULSIVE CORE FORCES. Carl Werntz, School 
of Physics, University of Minnesota, Minneapolis, 
Minnesota (Received August 3, 1960). 


A variational calculation of the binding energy 
of the triton has been carried out using the Gar- 
tenhaus potential. The results indicate that this 
potential leads to an unbound ground state of the 
three-nucleon system; this result is attributable 
to the even- parity tensor potential which is re- 
latively large in magnitude compared to the 
weakly attractive even-parity central potential. 
Since this property is also a characteristic of 
the Signell- Marshak potential, it too should lead 
to an unbound triton. 


COMPARISON OF INELASTIC SCATTERING 
FROM Sm?** WITH COULOMB EXCITATION 
THEORY. E. M. Bernstein and E. Z. Skurnik, 
Institute for Theoretical Physics, University of 
Copenhagen, Copenhagen, Denmark (Received 
September 9, 1960). 


In order to check the theory of electric quad- 
rupole Coulomb excitation, accurate measure- 
ments have been made of differential cross sec- 
tions for inelastic scattering from the first ex- 
cited state of Sm***. Protons, deuterons, and 
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alpha particles of energies in the region of 4 

Mev were used as the bombarding particles. In 
addition to the angular distribution of inelasti- 
cally scattered deuterons, measurements were 
also made with protons at two scattering angles 
and at a backward angle with deuterons and alpha 
particles of different incident energies. The data 
are found to be in excellent agreement with the 
semiclassical description of the Coulomb excita- 
tion process. 


DECAY OF ,,Er’® (3.1 hr). H. A. Grench and 
S. B. Burson, Argonne National Laboratory, 
Argonne, Illinois (Received July 8, 1960). 


Sources of ,,.Er'™ (3.1 hr) were produced by 
the (m,2n) and (y,m) reactions. A study of the 
gamma-ray spectrum by means of scintillation 
coincidence spectrometry indicated 32 gamma 
rays. All but three of these are fitted into a 
tentative decay scheme with levels at 0, 211, 
585, 826, 945, 1165, 1253, 1315, 1395, 1450, 
(1604), 1700, and 1830 kev. The data are con- 
sistent with a 3/2” spin assignment to the Er 
ground state and 7/2” and 1/2* assignments to 
the ground state and 211-kev state of Ho, respec- 
tively. The 826-kev level in Ho probably has 
5/2” character. 


Al?"(n, a)Na** CROSS SECTION AS A FUNCTION 
OF NEUTRON ENERGY. H. W. Schmitt and 
J. Halperin, Oak Ridge National Laboratory, Oak 


Ridge, Tennessee (Received September 16, 1960). 


The cross section for the Al?"(n,a)Na™ reaction 
has been measured as a function of neutron en- 
ergy in the range 6.1 <E, < 8.3 Mev and at 14.8 
Mev. Measurements were made relative to the 
fission cross section of U***; activation techniques 
were used to determine the number of Al”"(n, a) 
events. While a number of peaks and valleys ap- 
pear in the cross section versus energy curve, 
there is a general increase in cross section with 
increasing energy consistent with the Coulomb 
penetrability of the alpha particle. 


CHARGE -EXCHANGE SCATTERING OF NEGA- 
TIVE PIONS BY HYDROGEN AT 230, 260, 290, 
317, AND 371 Mev. John C. Caris, Robert W. 
Kenney, Victor Perez-Mendez, and Walton A. 


536 


Perkins, Il, * Lawrence Radiation Laboratory, 
University of California, Berkeley, California 
(Received August 29, 1960). 


The differential cross section for charge- 
exchange scattering of negative pions by hydro- 
gen has been observed at 230, 260, 290, 317, 
and 371 Mev. The reaction was observed by de- 
tecting one gamma ray from the 7° decay with a 
scintillation-counter telescope. A least-squares 
analysis was performed to fit the observations to 
the function 


do : 
do” D4. s00000 
l=1 
in the c.m. frame. The best fit to our experi- 


mental measurements requires only s- and 
p-wave scattering. The results (in mb) are: 





E(Mev) a, a, a, 

230+9 2.5040.10 1.3940.15 2.73+0.28 
260+7 2.02+0.08 1.754+0.14 2.15+0.22 
290+9 1.45+0.06 1.80+0.10 1.89+0.18 
31748 1.40+0.06 1.85+0.10 1.50+0.17 
371+9 1.084+0.05 1.6340.08 1.18+0.12 


The least-squares analysis indicates that d-wave 
scattering is not established in this energy range. 


*Now at Lawrence Radiation Laboratory, Livermore, 
California. 


LIFETIME OF THE NEUTRAL PION. V. Glaser, 
CERN, Geneva, Switzerland, and R. A. Ferrell, 
University of Maryland, College Park, Maryland 
and CERN, Geneva, Switzerland (Received Sep- 
tember 6, 1960). 


As Primakoff has noted, the phenomenological 
coupling constant of the neutral pion with the elec- 
tromagnetic field can be investigated by consider- 
ing the photoproduction of neutral pions in an ex- 
ternal Coulomb field. This is the inverse of the 
usual two-photon decay (one of the photons being 
provided by the external field). The relationship 
between the cross section and the free lifetime 
of the 7° is derived. Although the total cross sec- 
tion is small, it is found at high energy that the 
differential cross section is strongly peaked near 
the forward direction. The peak cross section is 
proportional to the fourth power of the photon en- 
ergy. It is this feature which makes possible an 
experimental determination of the lifetime by the 
photoproduction method to an accuracy of about 
ten percent. A minimum photon energy of one 














ir 
S 
1- 


e 











VoLUME 5, NUMBER 11 


PHYSICAL REVIEW LETTERS 


DECEMBER 1, 1960 





Gev is required to avoid uncertainties in the nu- 
clear form factor. A higher photon energy would 
be necessary only if the 7° mean life is greater 
than 5x10-'? sec. The backgrounds to be expected 
from nuclear photoproduction are estimated and 
found to be sufficiently small. In particular, the 
interference between the coherent nuclear 7° photo- 
production and the Primakoff process is not ex- 
cessive. 


RELATIVISTIC FORMULA FOR THE SPIN COR- 
RELATION COEFFICIENT Cxp. D. W. L. 
Sprung, Department of Mathematical Physics, 
University of Birmingham, Birmingham, England 
(Received September 27, 1960). 


A formula due to Stapp for the spin correlation 
coefficient in nucleon-nucleon scattering is cor- 
rected. Formulas for the other triple scattering 
and unpolarized incident beam correlation para- 
meters are included. The formulas are applicable 
when the nonrelativistic scattering matrix formal- 
ism of Wolfenstein is used to analyze triple-scat- 
tering experiments at higher energies where re- 
lativistic effects are not negligible. 

The definition of the correlation coefficients in 
the relativistic case is discussed. 


DETERMINATION OF PION-PION SCATTERING 
AMPLITUDES SATISFYING DISPERSION RELA- 
TIONS AND UNITARITY. John W. Moffat, RIAS, 
Baltimore, Maryland (Received September 2, 
1960). 


A method is developed for determining the 
partial-wave scattering amplitude in terms of 
the unitarity condition and the known branch cuts 
and poles of the inverse amplitude. The method 
is applied to the problem of pion-pion scattering 
and an implicit solution to the pion-pion partial- 
wave amplitude is derived for any angular mo- 
mentum state and for both elastic and inelastic 
scattering. With the aid of this solution the low- 
energy resonance behavior of the pion-pion 
scattering system is studied by neglecting all 
inelastic processes and concentrating on S and 
P waves. It is found that a P-wave resonance 
with a position and width required by nucleon 
electromagnetic structure can be determined in 
terms of two parameters. An iteration proce- 
dure is described that is applicable when the P 
wave dominates the equations and this procedure 


determines the contribution of the unphysical 
cut. The first iteration of the unphysical cut is 
numerically integrated on the IBM 709, and the 
results show that the shift of the resonance 
position due to the unphysical branch cut can be 
neglected. 


BOSON CURRENTS IN THE THEORY OF WEAK 
INTERACTIONS. W. B. Zeleny* and A. O. Barut, 
Department of Physics, Syracuse University, 
Syracuse, New York (Received May 24, 1960). 


A possible new method of introducing boson 
currents into the weak interaction is suggested. 
The Bose fields are assumed to obey the first 
order Kemmer equation. In this way, the boson 
and fermion currents enter the interaction in a 
symmetric manner. Three coupling constants 
are introduced according to the degree of iso- 
topic spin and strangeness symmetries of the 
currents. The interaction so obtained gives rise 
in a natural way to the |AT'| =3 selection rule 
and its violation, and also to fast pionic modes 
and slow leptonic modes of hyperon decay. First 
order calculations are performed for Kj.*, 1, 
and 7’ decays. Calculations of the decay rates 
yield sufficient information to determine the ap- 
proximate magnitudes of the effective coupling 
constants. The first order calculations are at 
variance with the experimentally uniform distri- 
bution of 7 events in the Dalitz plot. 


*Now at The University of Sydney, Sydney, N.S.W., 
Australia. 


STOCHASTIC DYNAMICS OF QUANTUM ME- 
CHANICAL SYSTEMS. E. C. G. Sudarshan, De- 
partment of Physics and Astronomy, University 
of Rochester, Rochester, New York, P. M. 
Mathews, Department of Physics, University 
of Madras, Madras, India, and Jayaseetha Rau, 
Department of Physics, Brandeis University, 
Waltham, Massachusetts (Received August 15, 
1960). 


The most general dynamical law for a quantum 
mechanical system with a finite number of levels 
is formulated. A fundamental role is played by 
the so-called “dynamical matrix” whose properties 
are stated in a sequence of theorems. A necessary 
and sufficient criterion for distinguishing dynami- 
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cal matrices corresponding to a Hamiltonian time 
dependence is formulated. The non-Hamiltonian 
case is discussed in detail and the application to 
paramagnetic relaxation is outlined. 


SPIN- MOMENTUM CORRELATIONS IN POSI- 
TRON-ELECTRON SCATTERING. C. Fronsdal 
and B. Jaksic, CERN, Geneva, Switzerland (Re- 
ceived September 27, 1960). 
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The imaginary part of the fourth order Bhabha 
(e+ - e-) scattering matrix element interferes with 
the (real) second order matrix element, to pro- 
duce a sixth order dependence of the cross sec- 
tion on the spin of one of the particles (after sum- 
ming over the spins of the other three particles). 
The process of extracting the imaginary part of 
the fourth order matrix element is presented in 
some detail for one of the graphs (vacuum polari- 
zation). 
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